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We review recent advances in the theory of the three-dimensional dilute homogeneous Bose gas 
at zero and finite temperature. Effective field theory methods are used to formulate a systematic 
perturbative framework that can be used to calculate the properties of the system at T = 0. 
T he p erturbative expansion of these properties is essentially an expansion in the gas parameter 
\/na 3 , where a is the s-wave scattering length and n is the number density. In particular, the 
leading quantum corrections to the ground state energy density, the condensate depletion, and 
long-wavelength collective excitations are rederived in an efficient and economical manner. We also 
discuss nonuniversal effects. These effects are higher-order corrections that depend on properties 
of the interatomic potential other than the scattering length, such as the effective range. We criti- 
cally examine various approaches to the dilute Bose gas in equilibrium at finite temperature. These 
include the Bogoliubov approximation, the Popov approximation, the Hartree-Fock-Bogoliubov 
approximation, the ^-derivable approach, optimized perturbation theory, and renormalization 
group techniques. Finally, we review recent calculations of the critical temperature of the di- 
lute Bose gas, which include l/Af-techniques, lattice simulations, self-consistent calculations, and 
variational perturbation theory. 
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I. INTRODUCTION 

The remarkable realization of Bose- Einstein conden- 
satio n (BEC) of trapped alkali atoms ll Anderson et all 
1199.4 iBradlev et flJ] . ll995UDavis et oIl . ll995l) has created 
an enormous interest in the properties of the weakly in- 
teracting Bose gas. Although the experiments are car- 
ried out in magnetic and optical harmonic traps, the 
ho mogeneous Bose gas ha s also received renewed inter- 
est l|Shi and GrimnlTl998(l . The homogeneous Bose gas 
is interesting in its own right and it may prove useful to 
go back to this somewhat simpler system to gain insight 
that carries over to the trapped case. 

Bose-Einstein condensation has a very long history 
dating back to the early days of qu antum mechanics and 
the p apers by Bose and Einstein ijBosel Il924t lEinsteiri 
1925). The canonical example of a system that exhibits 
BEC is liquid 4 He. At very low athmospheric pressure, 
4 He becomes superfluid below a temperature of 2.17K, 
which is called the A-point. The fluid then consists of 
two components, namely a normal and a superfluid com- 
ponent. The superfluid component has zero viscosity and 
it has the remarkable property that it can flow through 
narrow tubes without friction. The phase transition from 
a normal fluid to a superfluid is in modern terminology 
described as spontaneous breaking of the global U(l)- 
symmetry and the occurence of a condensate of atoms 
residing in the zero-momentum state. This description 
also explains the experimental fact that the spectrum is 
gapless and linear in the long- wavelength limit; the ex- 
istence of a gapless and linea r dispersion relation follows 
from the Goldstone theorem IjGoldstonel Il96ll) . It turns 
out that 4 He is strongly interacting and this reduces the 
density of the condensate to a rather small fraction of the 
total number density. At very low temperature and pres- 
sure, the condensate density is approximately 10% of the 
total number density. Thus the condensate density and 
superfluid density are very different at low temperature. 
Liquid 4 He consists of extended objects with very com- 
plicated nonlocal interactions. Moreover, the fact that it 
is strongly interacting makes it impossible to apply per- 
turbative methods. This has led to the search for weakly 
interacting Bose gases. The trapped alkali gases are such 
systems and the advantage of them over liquid 4 He is that 
they behave like systems of point particles with simple 
local interactions. These gases have therefore become 
a very active field of research in the past decade. In- 
tereste d readers may consult the books by Pethick and 
Smith llPethick and Smithl.l2002ft. and by Pitaevskii and 
Stringar Hjplt^ievskTi^uidSt^nTariL 120031) . as well a s the 
review pape rs by Da l fovo e t al. l|Dalfovo et all Il999j) and 
by Leggett tLeeeeti EfflOlll . 

The concept of Bose-Einstein condensation has been 
applied in many areas of physics other than 4 He and 
thus a thorough understanding of it is important. For 
instance, many properties of a superconductor can be 
understood in terms of a condensate of pairs of electrons 



with opposite spin and momenta. Similarly, many prop- 
erties of the QCD vacuum can be understood on the basis 
of a condensate of quark- antiquark pairs with zero total 
momentum. This condensate is called the chiral conden- 
sate and in massless QCD, it breaks the chiral symmetry 
of the QCD Lagrangian. The pions are interpreted as the 
corresponding Goldstone particles. 

The homogeneous Bose gas at zero temperature 
has been studied extensively for ove r 50 years startin g 
with the classic paper by Bogoliubov (Bo goliubovlll947(l . 
At zero temperature, the quantum loop expansio n es- 
sentially is an expansion in the gas parameter V na 3 , 
where n is the number density and a is the (positive) 
s-wave scattering length. The leading quantum correc- 
tions to the chemical potential, energy density and speed 
of sound, were calculated many years a go by Lee and 
Yang using the pseudo-potential method l)Lee and Yand . 
Il957j) . Part of the second quantum c orrec t ion to the 
energy density was o btained by Wu ifWu] Il959h . by 
Hugenholz a nd Pin es (Hugenhol z and Pinesl 11959). and 
by Sawada ijSawadal Il959|) . Only recently a com- 
plete two-l oop result has been obt ained by Braaten 
and Nieto ijBraaten and Nietd. Il997|) . The result de- 
pends, in addition to the scattering length a, also on 
an energy-independent term in the scattering amplitude 
for 3 — > 3 scattering. The fact that physical quanti- 
ties depend on properties other than the s-wave scatter- 
ing len gth was already poin t ed ou t by Hugenholz and 
Pines I Hugenholz and Pinesl Il959j) . These effects are 
called nonuniversal effects. In the past decades effective 
field theory has been established as an ideal tool for sys- 
tematically calculat i ng low-energy observables of a phys - 
ical system (|GeorgiL Il993t iKaolanL Il995t lLepagelll989T) . 
Effective field theory methods have also proven very use- 
ful for calculating higher order corrections in powers of 
V na 3 as well as nonuniversal effects i n the weakly inter- 
acting Bose gas l)Braaten et oil I2001|l 

Finite temperature corrections to the pressure of 
a dil ute Bose gas were first calculated by Lee and 
Yang IjLee and Yanel . 119581) . Performing an expansion of 
the pressure about zero temperature, they showed that 
the leading term goes as T 4 . This shows that the thermo- 
dynamics at low temperature is completely determined 
by the linear part of the spectrum. Similarly, the low- 
temperature expansio n of the numbe r dens ity was cal- 
culated by Glassgold ijGlassgold et a^J . ll960|) . and shows 
a T 2 behavior. These calculations were all based on the 
Bogoliubov approximation and are therefore valid only at 
low temperature, where the depletion of the condensate 
is small. The first approach that takes into account the 
excited states and thus can be used in the whole tem- 
perature range up to T r is the P opov theory ijPopovl 
I1983L Il987t IShi and Griffinl Il998t) . Very recently, an 
improve d many-body T-matrix approximation was de- 
veloped ijAl Khawaia et all l200l lAndersen et al\ . \200%\ 
which does not suffer from the infrared divergences in 
lower dimensions that plague the Bogoliubov and Popov 
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approximations. While the approach can be applied in 
any dimension, its main application is in one and two 
dimensions since the predictions in three dimensions are 
very similar to those of established approaches. For in- 
stance, it predicts the same zero-temperature depiction 
of the condensate as the Bogoliubov approximation and 
the same critical temperature as the Popov approxima- 
tion. 

The latest extensive review on t he homogeneous 
Bose gas was written by Shi and Griffin l)Shi and G riffin. 
1998) five years ago. In the mean time, significant 
progress has been made, and a new review paper 
summarizing recent advances is appropriate. There 
will of course be an overlap between the review by 
Shi and Griffin and the present paper, but some 
new material is covered. This i ncludes nonuniver- 
sal effects l|Braaten l20fllfl. renorma li zation 

group calculation s ll Andersen and Strickland! . Il999t 
iBiilsma and Stool Il996al iMetikas and Alberl 120021). 



some of the variational approaches IjLundh and Rammerl 
|2002() and calculations of the c ritical temperature 



T r (Arnold and Tomasik, 2000: Arnold and Moore, 


20011 
2000; 


lArnold et all l2002| iBavm et all 11999.2001 


Braatcn and Radescul l2002t iDavis and Morgan, 




Griiter et all Il997t iHp lzmann and Krauth. 1999 


Kashu 


rnikov et all 120011 iKastenind. l2003albl: iKleinert 


2003a 


iLedowski et a/1 120031 Ide Souza Cruz et all 12001 


2002). 


Traditionally, the theory of dilute Bose gases 



is presented using creation and annihilation operators 
as well as normal and anomalous self-energies. The 
present review is different in this respect since we are 
using modern functional methods and only occasionally 
we will present the material in terms of normal and 
anomalous averages. It also implies that we are working 
with Lagrangians and actions, rather than with Hamil- 
tonians. We are largely going to base o u r discussion on 
effective field theory methods llGeorgil 119931 iKaplanl 
Il995t iLepagel Il989|) . The virtues of effective field 
theory are that it allows one to systematically calculate 
physical quantities efficiently by taking advantage of the 
separation of length scales in a particular problem. We 
shall return to these issues in some detail later. 

Much of the research on the theoretical side on 
trapped Bose gases in recent years has been on dy- 
namical issues, e.g. condensate formation, damping of 
collective excitations, and collapse of the condensate. 
There is large body of literature on the noncquilibrium 
dynamics of the dilute Bose gas, but these topics de- 
serve a separate review. Similarly, following their ex- 
perimental realization, there has also been significant 
progre ss in the understanding of low-dimensiona l Bose 
gases ilAl Khawaia et all l2002t lAndersen et all l2002t 
iPetrov et all\200(t) . In order to limit the material, we 
will focus entirely on the three-dimensional Bose gas in 
this review. For the same reason, we also restrict our- 
selves to the spinless Bose gas. 

The paper is organized as follows. In Sec. II, we 



briefly review the ideal Bose gas at finite temperature. 
In Sec. Ill, we discuss the weakly interacting Bose gas 
at zero temperature. A perturbative framework using 
effective field theory methods is set up. This frame- 
work is used to rederive the leading quantum correc- 
tions to various quantities and discuss nonuniversal ef- 
fects. In Sec. IV, the weakly interacting Bose gas at finite 
temperature is reviewed. We discuss the Hartree-Fock- 
Bogoliubov approximation, the Bogoliubov approxima- 
tion, and the Popov approximation Then we consider the 
^-derivable approach and optimized perturbation theory, 
which are variational approaches. Finally, renormaliza- 
tion group techniques are reviewed. In Sec. V, we discuss 
the calculation of the critical temperature T c using a va- 
riety of different techniques. In Sec. VI, we summarize 
and conclude. Calculational details as well as notation 
and conventions are included in an appendix. 



II. THE IDEAL BOSE GAS 

In this section, we review the ideal Bose gas at fi- 
nite temperature. Although this is standard textbook 
material, we include a discussion to make the paper 
self-contained. In the remainder of the paper, we set 
Ti = 2m = ks = 1- Factors of h, 2m, and ks can be 
reinserted using dimensional analysis. 

Consider N bosons at temperature T in a box of vol- 
ume V, so that the number density is n — N/V. We im- 
pose periodic boundary conditions. We are always work- 
ing in the thermodynamic limit, meaning N, V — >■ oo in 
such a manner that n is fixed. 

At high temperature, where the thermal wavelength 
At = 2^/tF/T is much shorter than the interparticle spac- 
ing, the atoms behave classically and their statistics is not 
important. As the temperature is lowered, the atoms can 
be viewed as little wavepackets with extent of the order 
At- Bose-Einstein condensation takes place when the 
thermal wavelength of a particle is on the order of the 
interparticle spacing n -1 / 3 and the wavefunctions of the 
bosons start to overlap. The particles then accumulate in 
the zero-momentum state. At T = 0, all the particles re- 
side in this state. Since fermions behave very differently 
at low temperatures due to the Pauli exclusion principle, 
BEC is truly a quantum phenomenon. 

We can estimate Tj? by equating the thermal wave- 
length At with the average distance between the bosons 

n -l/3. 

2^M ~ n-^ . (1) 

The estimate for the critical temperature T® then be- 
comes 

T c ° ~ 47m 2 / 3 . (2) 
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Having estimated T c , we next discuss the thermo- 
dynamics in more detail. In the functional ap- 
proach to the imaginary-time formalism, the grand- 
cano nical partition function Z is given by a path inte- 
gral l)Negele and OrlandLll98^ : 



where the action S is given by 



S[ip*,i)] = J dr I d 3 x V*(x,t) 



(3) 



V(x,r) . (4) 



Here, (3 — 1/T and p is the chemical potential. The 
complex helds "0* and ip satisfy the standard bosonic pe- 
riodicity condition that 0(x, r) and , 0*( x j r ) are periodic 
in r with period /3. We next write the complex held in 
terms of two real helds a 



1 



ip = -j= (V>i + #2) 



(5) 



After inserting JSJ into the action Q, the corresponding 
Green's function or propagator is found to be 



D (u n ,p) 



1 



ojI + (p 2 - fi) 2 v p - v 



P - p w « 
,2 



(6) 



where tu n — 2nnT is the nth Matsubara frequency. The 
path integral © is Gaussian and can therefore be evalu- 
ated exactly. The result is 



_ - dr j d 3 x det £>~ 1 (w„ ,p) 



(7) 



where Dq 1 (uj n , p) is the inverse of the propagator JSJ and 
det denotes a determinant in the functional sense. The 
free energy density is given by 



TvlogD^^p) , 



(8) 
(9) 



where have used that log det A — Tr log A for any matrix 
A where Tr denotes the trace. Inverting, Eq. © and 
taking the trace, we obtain 



T = log [ujI + (p 2 - A*) 2 ] 



(10) 



In traditional approaches one uses the single complex field ip 
instead of two real fields. This is merely a matter of taste. See 
also the comments at the end of page 10. 



where the sum-integral means 



(11) 



The sum-integral involves a summation over Matsubara 
frequencies and a regularized integral over d = 3 — 2e 
dimensions. M is a renormalization scale that ensures 
that the integral has the canonical dimensions also for 
d 7^ 3. In the following, we absorb the factor M 2e in the 
measure and so it will not appear explicitly. We discuss 
the details further in the appendix. After summing over 
Matsubara frequencies, we obtain 



d d p 
{2n) d 



M)+Tlog 



(12) 



The hrst term inside the brackets is an infinite constant 
that is independent of temperature. It represents the 
zero-point fluctuations and can be removed by a vacuum 
energy counterterm AS. The second term is the standard 
finite temperature free energy of an ideal gas of nonrela- 
tivistic bosons. 

In the operator approach, the starting point would 
be the grand canonical Hamiltonian which corresponds 
to the action Q: 



d 3 p 
(27F 



(13) 



Using the Hamiltonian l|13l) to calculate physical quanti- 
ties, one does not encounter zero-point fluctuation terms 
since it has been normal ordered; The operator a p annihi- 
lates the vacuum. If it had not been normal ordered, one 
would find the same divergent terms as with the path- 
integral approach. This is the usual ambiguity of the 
quantization procedure in going from classical field the- 
ory to quantum field theory. 

We have replaced momentum sums by integrals over 
p. Due to the measure d d p, the integrand always vanishes 
at p = and the contribution from the ground state is 
not accounted for. If we denote the condensate density of 
particles in the lowest energy state by no and the number 
density of particles in the excited states by n ex , we have 
n = no + n GX . The density of excited particles is then 
given by minus the derivative of T with respect to the 
chemical potential: 



d d p 



n{p 2 - u) , 



where 
n(u>) - 



o/3w _ l 



(14) 



(15) 
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is the Bose-Einstein distribution function. The expres- 
sion H14JI makes sense only for fi < 0. If // > p 2 for some 
p* , the occupation number of the states with p < p* 
would become negative. Clearly, this is unphysical. Be- 
low the transition temperature, the number of particles 
in the excited states are given by the integral Ijl4|) with 
H = 0. One finds 



c(l) 

(4tt) 3 / 2 



T 3/2 



(16) 



where £(a;) is the Riemann zeta function with argument 
x. The critical temperature is the temperature at which 
all the particles can be accomodated in the excited states, 
that is n = n cx . This yields 



47T 



c(l) 



2/3 



(17) 



The quantity is called the degeneracy parameter. For 
an ideal Bose gas, the critical number density n c satisfies 
n c A^ = C(|)- We see that the estimate (|2J) is correct 

within a factor of [( (|)j 2 ^ 3 0.527 of the exact re- 
sult (XZt - Using Eqs. Ijltil) and l|17l) . the condensate den- 
sity no as a function of the temperature T can be written 
as 



n 



3/2' 



(18) 



The exponent 3/2 in Eq. 1|18|) is determined solely by 
the density of states. For an ideal Bose gas in a three- 
di mensional isotropic harmonic trap, the exponent is 
3 l|Dalfovo et oJ] . ll999h . 



III. WEAKLY INTERACTING BOSE GAS AT ZERO 
TEMPERATURE 



In this section, we discuss in some detail the weakly 
interacting Bose gas at T = 0. We begin with a de- 
scription using effective field theory and formulate a per- 
turbative framework that can be used for practical cal- 
culations. We then calculate the leading corrections in 
the low-density expansion to the energy density, deple- 
tion, and long- wavelength excitations. Finally, we discuss 
nonuniversal effects. 



iKaplanl Il995l iLepagel Il989|) . EFT takes advantage of 
the separation of scales in a system to make model- 
independent predictions at low energy. The effective 
Lagrangian C e g that describes the low-energy physics is 
written in terms of only the long-wavelength degrees of 
freedom and the operators that appear are determined by 
these degrees of freedom and the symmetries present at 
low energy. The effective Lagrangian generally includes 
an infinite tower of nonrenormalizable interactions, but 
they can be organized according to their importance at 
low energy; to a certain order in a low-energy expan- 
sion, only a finite number of operators contribute to a 
physical quantity and one can carry out renormalization 
in the standard way order by order in this expansion. b 
Since the coefficients of these operators encode the short- 
wavelength physics, we do not need to make any detailed 
assumptions about the high-energy dynamics to make 
predictions at low energy. 

In some cases, one can determine the coefficients of 
the low-energy theory as functions of the coupling con- 
stants in the underlying theory by a perturbative match- 
ing procedure. One calculates physical observables at low 
energies perturbatively and demand they be the same in 
the full and in the effective theory. The coefficients of the 
effective theory then encode the short-distance physics. 
If one cannot determine these short-distance coefficients 
by matching, they can be taken as phenomenological pa- 
rameters that are determined by experiment. 

A classic example of an effective field theory is chiral 
pertur bation theory which is a low-e nergy field theory for 
pions (jGasser and Leutwvlerl Il984[) . Pions are interact- 
ing particles whose fundamental description is provided 
by QCD. However, QCD is strongly interacting and con- 
fining at low energies and so perturbative calculations 
using the QCD Lagrangian are hopeless. So instead of 
using the quark and gluon degrees of freedom in QCD, 
one writes down the most general Lagrangian for the pi- 
ons, which are the relevant low-energy degrees of free- 
dom. The terms that appear in the chiral Lagrangian are 
determined by the global symmetries of QCD. The coef- 
ficients of the chiral Lagrangian cannot be determined as 
functions of the couplings and masses in QCD using per- 
turbative methods. They can in principle be determined 
using a nonperturbative method such as lattice gauge 
theory, but in practice they are normally determined by 
experiment. 

Nonrelativistic QED ijCaswell and Lepaeel Q.986) 
(NRQED) is an example of an effective field theory whose 
coefficients are tuned so they reproduce a set of low- 



A. Effective Field Theory 

Effective field theory (EFT) is a general approach 
that can be used to analyze the low-ener gy behavior of 
a physical system in a systematic way (|Georgl Il993l 



b The counterterms that are used to cancel the divergences in 
the calculations of one physical quantity are the same as those 
required in the calculations of another. Thus, having determined 
the counterterms once and for all, the effective theory can be used 
to make predictions about other physical quantities. 
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energy scattering amplitudes of full QED. NRQED is tai- 
lored to perform low-energy (bound state) calculations, 
where one takes advantage of the nonrelativistic nature of 
the bound states by isolating the contributions from the 
relativistic momentum scales. These are encoded in the 
coefficients of the various local operators in the effective 
Lagrangian. Traditional approaches involving the Bcthc- 
Salpeter equation do not take advantage of the separation 
of scales in bound state problems and are therefore much 
more difficult to solve. 

Landau's quasiparticle model for 4 He can also be 
viewed as an effective theory. In order to explain that the 
specific heat varies as T 3 for temperatures much smaller 
than T c , he suggested that the low- lying excitations are 
phonons with a linear dispersion relation. More gen- 
erally, he proposed a spectrum that is linear for small 
wavevectors and has a local minimum around p = po. 
This part of the spectrum behaves like 



tral part of the interaction is the magnetic dipole-dipole 
interaction. 

The central part of the potential consists of a short- 
range part with range xq and a long-range van der Waals 
tail. The latter goes as 1/x 6 as x — > oo. A typical two- 
body potential is shown in Fig. ^ A model potential of 
this kind is the sum of a hard-core potential with range 
xq and a van der Waals potential: 



V c (x) = 



-oo , X < Xq 

, x > x 



(21) 



where C$ and xq are constants. Another example is the 
hard-core square-well potential: 



V c (x) = 



1-00 , X < X c 
-Vq , Xq < X < X c 

, x < x c 



(22) 



2toq 



(19) 



where A, po, and mo are phenomenological parameters. 
Excitations near pg are referred to as rotons. Assuming 
that the elementary excitations are noninteracting, one 
can use spectrum to calculate the specific heat. Landau 
determined the parameters by fitting the calcu lated spe- 
cific heat to experimental data ijLandaul Il947|) . 

The weakly-interacting Bose gas is a system where 
effective field theory methods can be applied successfully. 
The starting point is the action: 



SWA] 




d a x tp*(x,t) 



J n2 

*m + V +M 
ot 



xip(x,t)-- I d d x / A'f(x,f)f(x',i) 



xF (x-x')V'(x,t)V'(x / ,<) + 



(20) 



Here, i/>*(x, t) is a complex field operator that creates 
a boson at the position x at time t, \i is the chemical 
potential, and Vb(x) is the two-body potential. The el- 
lipses indicate terms that describe possible interactions 
between three or more bosons. The chemical potential fi 
must be adjusted to get the correct number density n. 

The interatomic potential can be divided into a cen- 
tral part Vq(x) and a remainder. The central of the po- 
tential only depends on the separation x of the atoms and 
their electronic spins. It conserves separately the total 
orbital angular momentum and the total electronic spin 
of the atoms. The noncentral part of the interaction con- 
serves the total angular momentum, but not separately 
the orbital angular momentum and the total electronic 
spin of the atoms. An example of a term in the noncen- 



where xq and constants. This potential sustains 

a number of two-body bound states depending on the 
values of xq and x c . Many real potentials used in ex- 
periments sustain bound states and the ground state of 
the system is no longer a homogeneoues gas, but rather 
a state of clusters of atoms. However, if the scattering 
length is positive, the homogeneous Bose gas can exist as 
a long-lived metastable state. 

In Fig. |21 we have shown the Fourier transform V(k) 
of a typical short-range two-body potential with range 
xq. Since a true interatomic potential vanishes for large 
momenta k, one will never face ultraviolet divergences 
using it in actual (perturbative) calculations. 

In this paper we restrict our calculations to the spin- 
less Bose gas. In experiments with trapped alkali atoms 
(e.g. 7 Li, 23 Na 85 Rb, 87 Rb, and 133 Cs), the situation is 
more complicated. In all these atoms, there is a single s- 
electron outside closed shells. Consequently, these atoms 
have electron spin S — 1/2. The total spin of a colliding 
pair of alkali atoms is therefore either S = or S = 1. 
The central part of the potential V c (x) depends on the 
total spin and one refers to these potentials as the singlet 
and triplet potentials, respectively. The singlet potential 
is generally much deeper and sustains many more bound 
states than does the triplet potential. For example, the 
singlet potential of 87 Rb is deeper than the triplet poten- 
tial by more than one order of magnitude. The scattering 
lengths are denoted by a s and at- In the case of atomic 
hydrogen, they have been calculated from first principles; 
a s = 0.3cto and a t — 1.3ao, where ao is the Bohr radius. 
Similarly, the coefficient of the van der Waals tail h as 
also been determined; Cq = 6.499ao l|Yan et all 119961. 

One can associate a natural length scale I with any 
atomic potential V c (x). For a short-range potential, this 
length is the range x itself. For a long-range potential, 
it is a little more complicated. The length scale asso- 



ciated with the van der Waals tail is I 



vdW 



C, 



1/4 



111 
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this case, the length scale I is either the range or Z v dw> 
whichever is larger. For a generic potential, the low- 
energy observables such as the scattering length and the 
effective range are of the order I. There is nothing that 
forbids the these quantities to be much larger than I, but 
it is unnatural and typically it requires fine-tuning of one 
or more parameters in the potential. In the case of the 
alkali atoms, / v dw is much larger than the range of the 
short-range part of the potential and is thus the natural 
length scale for these atoms. The spin-singlet scatter- 
ing length for 85 Rb is a s = +2800ao and is more than 
one order of magnitude larger than I = Z v dw = 164ao- 
This can be viewed as a fine-tuning of the mass of the 
atom. This can be seen from the fact that the mass of 
87 Rb is only 2.3% larger and the spin-singlet scattering 
length has the more natural value a s = +90.4ao. One 
can also obtain unnaturally large scattering lengths by 
tuning an external parameter in experiments. One way 
of doing this, is to t une an external magnetic field to a 
Feshbach resonance llFeshbachllT96l . This is currently 
receiving a lot of attention both theore tically |Stwallevt 
ll976tlTiesinga. Verhaar and Stooi fi993) and experimen- 
tally (ftnouve et a/.LIl998(l . 

In the remainder of this work, we do not specify the 
magntiude of the scattering length. We only demand the 
diluteness condition be satisfied, namely that 



V(x) 



ncr < 1 



(23) 



We now return to the action (|20|l . which is invariant un- 
der a global phase transformation 



*V(x,i) 



(24) 



The global U(l) symmetry reflects the conservation of 
atoms. It also ensures that the number density n and 
current density j satisfy the continuity equation 




FIG. 1 Typical behavior of a two-body potential Vq(x). 




Vj = 



(25) 



FIG. 2 Solid line: Typical behavior of the potential in mo- 
mentum space V(k). Dotted line: Fourier transform of a 
contact potential. 



where the dot denotes differentiation with respect to 
time. 

The nonlocal evolution equation that follows from 
the action l|2U|) . cannot always be used for practical cal- 
culations, but can be replaced by a local one. Suppose 
we are interested in the properties of the system at mo- 
menta k such that the de Broglie wavelength 1/fc is much 
longer than the range of the interatomic potential Vb(x). 
The interactions therefore appear pointlike on the scale 
of the de Broglie wavelength, and they can be mimicked 
by local interactions The parameters of these local inter- 
actions must be tuned so that they reproduce low-energy 
observables to sufficient accuracy. However, if the poten- 
tial is long range, the scattering amplitude depends in a 
nonanalytic way on the wave vector k characterizing the 
incoming atoms in the CM frame. Such behavior can- 
not be reproduced by local operators. If Vq(x) falls of 



like 1/x 6 , this nonanalytic beha vior enters first at order 
k 4 ijBraaten and Hammeit 12003^ . 

The effective Lagrangian for the bosons can be con- 
structed using the methods of effective field theory. Once 
the symmetries have been identified, one writes down the 
most general local effective Lagrangian consistent with 
these symmetries. At zero temperature, the symmetries 
are Galilean invariance, time-reversal symmetry, and the 
global phase symmetry l|24|) . These symmetries severely 
restri ct the possible terms in t he effective action. One 
finds l)Braaten and Nietd. [l997^ : 



S[p,il>] 



dt / d d x { tp* 



■ 9 w2 
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-\g {r^f - \h [v (rn 2 - § (r^f + 



(26) 



where g, h, and gs are coupling constants that can be de- 
termined by a matching procedure. The dots denote op- 
erators that are higher order in the field ip or its derivative 
V^> and respect the symmetries. One demands that the 
effective field theory represented by the action (|2l)|) re- 
produces a set of low-energy observables to some desired 
accuracy. An example is the coefficients of the expansions 
in V na 3 of the ground state energy density. Another ex- 
ample coefficients in the low-energy expansions for the 
scattering amplitudes in the n-body sector. When the 
coefficients in the action <|26|) have been determined to 
some accuracy in a low-energy expansion, effective field 
theory guarantees that all other observables can be de- 
termined with the same accuracy. 

The quantum field theory defined by the action 126|) 
has ultraviolet divergences that must be removed by 
renormalization of the parameters /z, g, h, g^... . They 
arise because we are treating the interactions between 
the atoms as pointlike down to arbitrarily short dis- 
tances. For instance the operator g(ij}^ip) 2 can be thought 
of as a contact potential with strength g; Vo(x — x') = 
gS(x — x'). The Fourier transform is then a constant in 
momentum space; Vb(k) = g. This is illustrated is Fig. [3 
where the dotted line shows U(k). Thus U(k) does not 
vanish for large momenta and this is the reason why one 
encounters ultraviolet divergences in the calculation of 
Feynman diagrams. In order to make the theory well 
defined, we must introduce an ultraviolet cutoff. This 
is indicated in Fig. [21 where we exclude wave numbers 
k > M in momentum integrals. 

If we use a simple momentum cutoff M to cut off 
the ultraviolet divergences, there will be terms that are 
proportional to M n , where n is a positive integer. There 
are also terms that are proportional to log(Af). The co- 
efficients of the power divergences depend on the method 
we use to regulate the integrals, while the coefficients of 
log(M) do not. Thus the power divergence are artifacts of 
the regulator, while the logarithmic divergences represent 
real physics. I n this paper, we will be using dimensional 
regularization i j't Hooft and Veltmanl . Il972f) to regulate 
infrared as well as ultraviolet divergences in the loop in- 
tegrals. In dimensional regularization, one calculates the 
loop integrals in d = 3 — 2e dimensions for values of e 
for which the integrals converge. One finally analytically 
continues back to d = 3 dimensions. In dimensional regu- 
larization, an arbitrary momentum scale M is introduced 
to ensure that loop integrals have their canonical dimen- 
sions also away from three dimensions. This scale can be 
identified with the simple momentum cutoff mentioned 
above. Advantages of dimensional regularization are that 
it respects symmetries such as rotational symmetry and 
gauge invariance. Dimensional regularization sets power 



divergences to zero and logarithmic divergences show up 
as poles in e. With dimensional regularization, the only 
ultraviolet divergences that require explicit renormaliza- 
tion are therefore logarithmic divergences. This simpli- 
fies calculations significantly as we shall see. In fact all 
the divergences encountered in the one-loop calculations 
that we present here, can be removed by the renormal- 
ization of /x, g, and the vacuum energy £ and in three 
dimensions, these are power divergences. Thus no ex- 
plicit renormalization is required. 

We now return to the determination of the parame- 
ters g,h,gs... . We follow ijBraaten et all l2~00l[) and de- 
termine them by demanding that the effective field the- 
ory l|26|) reproduces the low-momentum expansions for 
the scattering amplitudes in the vacuum for 2^2 scat- 
tering, 3 — > 3 scattering, etc. to some desired accuracy. 

To calculate the coupling constant g, consider the 
scattering of two atoms in the vacuum with initial wave 
numbers ki and k2, and final wave numbers k^ and k 2 . 
The probability amplitude for the 2 — > 2 scattering is 
given by the T-matrix element T. The tree-level contri- 
bution to T comes from the leftmost diagram in Fig. [3] 
and is given by 



-2<7, 



(27) 



where the subscript indicates the number of loops. The 
quantum corrections to the tree-level result come from 
the loop diagrams in Fig. |3| The leading quantum cor- 
rection comes from the one-loop diagram and reads 



. . „ 2 f duo f d d k 1 

T x {q) = -2* 5 J -J {27r)dLu _ k , 



is 



(W - k\ - k\) + (k - ki - k 2 ) 2 + IE 



(28) 



where q = ^|ki — k&l, and we have used that the free 
propagator A.q(u>, k) in vacuum corresponding to the ac- 
tion (123) is 



A (w,k) 



-fc 2 



(29) 



The integral over to is performed using contour integra- 
tion. After changing variables k — * k + |ki -(- ki2 1 /2, we 
obtain 



d d k 



1 



{2TT) d k 2 - q 2 



ie 



(30) 



Note that the integral over k is linearly divergent in the 
ultraviolet. This divergence is set to zero in dimensional 
regularization. Using dimensional regularization, the re- 
sult of the integration of k can be written as 



Tx{q) 



g 2 M 2t V 



(47r) d / 2 



(31) 
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where T(x) is the gamma function. The limit d 
regular and Eq. (|31|1 reduces to 



: g 2 q 

'167T 



(32) 



This expression is simply Fermi's golden rule. The quan- 
tum corrections from higher orders are given by the di- 
agrams like the two-loop graph in Fig. [J] They form a 
geometric series that can be summed up exactly and the 
exact 2 — > 2 amplitude becomes 



2g 2 



9 + TM 



(33) 



The scattering amplitude for 2 — > 2 scattering in the 
underlying theory described by the action l|20[l can be 
calculated by solving the two-body scattering problem in 
the potential Vn(x — x'). The contribut ion from s-wave 
scattering is l)Landau an d Lifshitz, 1980) 



T 



16tt 



lSa(q) sin[<5 (<?)] 



(34) 



where So(q) is the s-wave phase shift. We next write the 
low-momentum expansion as follows 



be shown that the ladder diagrams are of equal impor- 
tance and must be summed. The summation of these 
diagrams can be expressed in terms of an effective in- 
teraction r that sati sfies an integra l equation which also 
invol ves Vo(x — x') ()Beliaevl Il958l iFetter and Waleckal 
Il97l|) . The interatomic potential appearing in T can be 
eliminated in favor of the scattering amplitude for two- 
particle scattering. In the low- momentum limit, the ef- 
fective interaction reduces to the 8™. The mean-field 
self-energies and therefore the first order propagator and 
spectrum can be expressed in terms of the effective inter- 
action. Thus the spectrum reduces in the low-momentum 
limit to the Bogoliubov spectrum. 



B. Perturbative Framework 

We next discuss the perturbative framewo rk set up 
by Braaten and Nieto ijBraaten and NietoLll997() that can 
be used to systematically calculate the low-energy prop- 
erties of a weakly interacting Bose gas. 

We first parameterize the quantum field in terms 
of a time-independent condensate v and a quantum fluc- 
tuating field 0; 



= v + . 



(37) 



qcot[S a (q)] = 



1 1 2 



(35) 



The fluctuating field can be conveniently written in 
terms of two real fields: 



This expansion defines the scattering length a and the 
effective range r s . Using the identity e lx sin(x) = 
l/(cot(x) — i), we can expand the T-matrix in powers of 
momentum q. Matching the expressions l|33l) and l|34|) 
through first order in the external momentum q us- 
ing i|35|). we obtain 



g = 87ra 



(36) 



The parameter h can be determined by going to the next 
order in the low-momentum expansion. We will do this 
in subsec. Illl.FI Similarly, the parameter 53 can be de- 
termined by solving the three-body scattering problem in 
the potential Vb(x— x') and demand that the 3 — > 3 scat- 
tering amplitudes in the full and in the effective theory 
be the same at low momentum. 






FIG. 3 Diagrams contributing to the 2 — » 2 scattering ampli- 
tude. 



i> = -j= C0i +#2) 



(38) 



Substituting Eqs. JSZJl and (jSHJ into Eq. iffify. the action 
can be decomposed into three terms 

S[v,i/>i,ip2] = S[v] + S ilee [v,ipi,ip 2 ] + SintfajV'ijV'a] ,(39) 

where we have indicated that the action depends on v as 
well, and switched to the variables ip\ and ip2 instead of 
ijj* and -0- S[v] is the classical action 



S[v] 



dt I d d x 



1 



[IV 



-gv 



while the free part of the action is 



<Sfreeb,V>l! "02 



dt / d d x 



^ (01-02 - 0i0; 



+ ^01 (V 2 + X) 0! + ^0 2 (V 2 + Y) 02 



where 



(40) 



(41) 



Traditionally, the starting point has been the inter- 
atomic potential Vq(x — x'). At low densities, it can 



X 
Y 



fx-3gv , 
H - gv 2 . 



(42) 
(43) 
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The terms 3gv 2 and gv 2 in X and Y are often referred 
to as mean-field self-energies. The interaction part of the 
action is 



dt / (fx 



The sources in Eq. (|44|l are 

J = [m - 9 v2 ] v , 

Z = -gv . 



,(44) 



(45) 
(46) 



The propagator that corresponds to the free action 
S{ Tee [v,Tpi,ip 2 ] in Eq. (gTJ is 



uj 2 — e 2 (p) + ie 



p — Y ~iuj 
iuj p 2 — X 



(47) 



Here, p = |p|, where p is the wavevector, u is the fre- 
quency, and e(p) is the dispersion relation: 



e(p) = V(P 2 ~ X)(P 2 ~ Y) . 



(48) 



Note that one can diagonalize the matrix (|4T|) by a field 
redefinition, which is equivalent to a Bogoliubov transfor- 
mation in the operator approach. Such a field redefinition 
makes, however, the interaction terms much more com- 
plicated and increases the number of diagrams that one 
needs to evaluate. For practical purposes, we therefore 
stick to the above propagator. 

The partition function Z can be expressed as 
a path integral over the quantum fields ip\ and 
•02 l|Negele and OrlandLll988(l : 



(49) 



All the thermodynamic observables can be derived from 
the partition function Z. For instance, the free energy 
density T is given by 



T{p) = i 



V 



(50) 



where V is the spacetime volume of the system. The 
pressure V(fi) is 



V(fx) = -T{n) 



(51) 



The number density n is given by the expectation value 
(0*0) in the ground state: 



i(fj) = / D0i£>0 2 (0*VO e lS[v '^'^ ] 



(52) 



It can therefore be expressed as 
dF(n) 



i(p) = 



dp. 



(53) 



The energy density £ is given by the Legendre transform 
of the free energy density T: 



(54) 



The free energy T is given by all connected vacuum 
diagrams which are Feynman diagrams with no exter- 
nal legs. The sum of the vacuum graphs is indepen- 
dent of the condensate v. At this point it is convenient 
to introduce the thermodynamic potential fl(p,,v). The 
thermodynamic potential is given by all one-particle irre- 
ducible vacuum diagrams and can be expanded in number 
of loops: 

Si(p,v) = SloOhv) +fii(A*.") + n 2 (At,u) + ... , (55) 

where the subscript n denotes the contribution from the 
nth order in the loop expansion. If f2 is evaluated at a 
value of the condensate that satisfies the condition 



V = (0) , 



(56) 



all one-particle reducible diagrams (those that are discon- 
nected by cutting a single propagator line) vanish. Thus 
evaluating the thermodynamic potential at the value of 
the condensate that satisfies 1)56)1 . one obtains the free 
energy: 



T(pi) = ft (fx, v) + f2i(/x,w) + 2 (/i, v) + 



(57) 



Using Eq. the condition (f^TTf) reduces to (^2) = 

and (0i ) = 0. The first condition can be automatically 
satisfied by a suitable choice of the phase of ifi. The 
second condition is then equivalent to 



9fi(/x, v) 
dv 







(58) 



The value of the condensate that satisfies l|58|l is denoted 
by v. The free energy can also be expanded in powers of 
quantum corrections around the mean-field result T${p): 



= ^o(m) + ^i(m) + Mf*) + 



(59) 



The loop expansion l|57|l of does not coincide with 
the expansion (|59|l of F(p) in powers of quantum correc- 
tions because of its dependence on v. To obtain the ex- 
pansion of T in powers of quantum corrections, we must 
expand the condensate v in powers of quantum correc- 
tions: 



V — V + Vl + v 2 + ... , 



(60) 
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where vq is the classical minimum, which satisfies 



dv 



= 



(61) 



By expanding Eq. (|58|l about v , one obtains the quan- 
tum corrections v\, V2, ■■■ to the condensate. For instance, 
the first quantum correction v% to the classical minimum 
v is 



vi = 



dv 



d 2 Q (fi,v) 



dv 2 



The mean-field free energy density is 



(62) 



(63) 



Inserting l|6U|) into (|57|) and expanding in powers of 
v\,v 2 ..., we obtain the quantum expansion of the free 
energy density. The first quantum correction to the free 
energy density is 



^"i(/x) = Qi(/i,u ) , 



(64) 



the dispersion relation is approximately e(p) = p 2 + p, 
where the second term represents the mean-field energy 
due to the interaction with the condensed particles. 

We next comment on how traditional approaches fit 
into the more general perturbative framework discussed 
in this section. The starting point is the second quantized 
grand canonical Hamiltonian that includes a two-body 
potential Vo(x): 

H = J d d x ip*(x) [-V 2 - fi] V»(x) 

+i J d d x J d d x' ?/»*(x)V'*(x')Fo(x-x')V(x)V'(x') . 

(68) 

This Hamiltonian is nonlocal and is often replaced by a 
local one. This is done by approximating the true two- 
body potential by a local two-body interaction, whose 
strength g is tuned to reproduce the scattering length of 
Fo(x-x'). This yields 



U = / d d x \ r [- V 2 - /i] ij + X -g (</>* V) 2 \ ■ i (>') ! 



and the second quantum correction to the free energy 
density is 



T 2 (\x) = n 2 {fJ,,v ) + vi 



dn^p.v) 



dv 



1 2 d 2 n a {n,v) 



dv 2 



(65) 



Expressions for higher order corrections to the free energy 
can be derived in the same way. 

The value of the condensate v that minimizes the 
classical action l|40f) is given by vq — \J~p]~9- The linear 
term in Eq. I|44|l then vanishes since J = 0. At the mini- 
mum of the classical action, both the propagator l|47[> and 
the dispersion relation l|48|l simplify significantly since we 
also have Y = 0. At the minimum, these equations re- 
duce to 



D(u,p) = 



uj 2 - e 2 (p) +is \iu e 2 (p)/p 



p ~IU) 



(66) 
(67) 



The spectrum <|6 7 ll was first derived by Bogoliubov in 
1947 l|Bogoliubovl 11947) . The dispersion relation is gap- 
less and linear for small wave vectors. This reflects 
the spontaneous breakdown of the Z7(l)-symmetry. The 
dispersion relation changes from being linear to being 
quadratic in the vicinity of p = ^/2jl, which is called the 
inverse coherence length. For very large wave vectors, 



The grand canonical Hamiltonian is expressed in terms 
of creation and annihilation operators that satisfy the 
standard equal-time commutation relations. Bogoli- 
ubov's idea was to treat the p = momentum sepa- 
rately ijBogoliubovl 1 1 947}) . Since this state is macroscop- 
ically occupied, the creation and annihilation operators 
for bosons with p = commute to very good approx- 
imation. Thus they can be treated classically and be 
replaced by a constant which is the condensate density. 
This step is equivalent to splitting the quantum field ip 
into a condensate v and a fluctuating field ip, Eq. (|37|l . 

I n the Bogoliubov approximation ijBogoliubovt 
11947ft . one makes a quadratic approximation to the 
Hamiltonian by neglecting terms with three and four op- 
erators. Since the resulting Hamiltonian contains prod- 
ucts of two annihilation and products of two creation 
operators, it must be diagonalized by a canonical trans- 
formation. The resulting quasi-particle spectrum is then 
given by Eq. 1)6 7|) . 

In the Beliaev approximation ijBeliaevl Il958ft . one 
goes one step further by calculating the leading quantum 
corrections to the quasi-particle spectrum (|67|) . This is 
done by including all one-loop diagrams, that is all dia- 
grams up to second order in the interaction, in the self- 
energies and calculating the poles of the propagator. The 
correction to the Bo goliuboy spect rum l|6T|l was first cal- 
culated by Beliaev ( Beliaevt Il958ft and coincides with a 
leading-order calculation in the approach outlined here. 
We shall return to that calculation in Sec. 1111. Dl 

The perturbative framework has been formulated in 
terms of two real fields ipi and ip 2 - As noted in the 
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introduction, one has traditionally presented the the- 
ory in terms of normal and anomalous Green's functions 
Gn(uj,p) and Gi2(w,p) and the self-energies Sn(w,p) 
and Ei2(w,p). This formulation corresponds to using the 
fields ip* an d ip instead of tpi and -02- In cither formula- 
tion, one ends up with a 2 x 2 matrix for the propagator 
and self-energies and the amount of work to calculate 
most quantities is comparable. For readers who want 
to translate intermediate results to a more familiar lan- 
guage, we note that 



Sn(w,p) 



[Uu(uj,p) - n 22 (w,p)] 



(70) 
(71) 



where Hij(ui,p) are the components of the 2x2 self- 
energy matrix. 

We next comment on the Hugenholz-Pines theo- 
rem ijHugenholz and Pinesl Il95^) . The Hugenholz-Pines 
theorem ensures that the spectrum does not exhibit a 
gap and it has been proven to all orders in perturbation 
theory. It is simply the Goldstone theorem for a dilute 
Bose gas with a spontaneously broken continuous sym- 
metry. The Hugenholz-Pines theorem is normally given 
in terms of normal and anomalous self-energies: 

H = En(0,0)-Sia(0,0) . 

In terms of the self-energies Hij(cj,p), the theorem takes 
a particular simple form: 



n = n 22 (o,o) 



(72) 



C. Ground State Energy Density and Condensate Depletion 

In this subsection, we calculate the leading quan- 
tum correction to the ground state energy and the deple- 
tion of the condens ate. While these results are standard 
textbook material ijFetter and Waleckal Il97lh , it is in- 
structive to see how they are derived within the present 
framework. 

The mean-field thermodynamic potential f2o(Mj v ) is 
given by the terms in the classical action (|4()(l : 



The mean-field number density follows from differentiat- 
ing l)74|) with respect to u. The chemical potential in the 
mean-field approximation is then obtained by inversion: 

= gn ■ (75) 

The mean-field energy £ is easily from <|54[1 found to be 

= Aiian 2 . (76) 

The mean-field result for the energy dens ity was first ob- 
tained by Bogoliubov ijBogoliubovl Il947j) . 

The one-loop contribution to the thermodynamic 
potential f2(/x, v) is given by 



Qx((i,v) = i 



(77) 



where Zq is the path integral involving the quadratic 
quantum fluctuations around the mean field. 



_ ifdtf d 3 KdctD _1 (w,p) 



(78) 



where D^ 1 ^^) is the inverse of the propagator (|4T|> and 
det denotes a determinant in the functional sense. Using 
the fact that TrlogA = log det A for any matrix A, we 
obtain 



~ / n 1 f du> f d d k , , -i . 

n^v) = --^-y^iogdet^-w) 

+A 1 n . 



(79) 



Here we have added Aifi, which is the one-loop countert- 
erm. The counterterm is added to cancel the ultraviolet 
divergences that one encounters when evaluating the in- 
tegral in J75J). After integrating over uj using (|A. 1|) . we 
obtain 



QlM = \ J 70V(fc 2 -*)(fc 2 -iO 



(80) 



flo(fjt,v) = -fiv 2 + -gv 4 



(73) 



The minimum of the mean-field thermodynamic poten- 
tial is given by vq = \ffjijg and the mean-field free energy 
J-q is obtained by evaluating Eq. I|73|) at the minimum: 



25 



(74) 



The one-loop contribution T\ to the free energy is ob- 
tained by evaluating Eq. Ij80|l at the classical minimum, 
where Y = and X = —2/j,. The one- loop free energy 
then reduces to 



' 2g ' 2 J (2?r) 
-^ + i/ ,-i(2/i) + A 1 ^), 



(81) 
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where the integral J mj „(A) is defined in the appendix and 
AiJF is the one-loop counterterm. The integral /o,-i(A) 
has quintic, cubic, and linear ultraviolet divergences that 
are set to zero in dimensional regularization. The coun- 
terterm A\T is therefore zero and the limit d — > 3 is 
regular. We obtain 



2g 



4VW 
15tt 2 



(82) 



It might be useful to see how the renormalization proce- 
dure works with a simple ultraviolet cutoff. In that case, 
the one-loop contribution to the free energy T\ can be 
written as 



^i(At) = 



M 



d 3 k 
(27T) 5 



k^k 2 + 2/i + AiJT , 



(83) 



where the integral is calculated in d = 3 dimensions and 
the superscript M indicates that |k| < M has been im- 
posed. We can now rewrite this as 



footnote |b| The renormalized one-loop free energy can 
then be written as 



•^o+i (m) 



1 f M d 3 k r 2 
2- 9 + 2j WT^ k) - k 

,2 



2/c 2 



(89) 



oo. Evaluating the 



We can now take the limit M 
integral, we recover Eq. JS^J. 

Finall y, we mention that the pseudo-potential 
method ((Huang and YaneL Il957|) is an alternative way 
of treating the ultraviolet divergences. One replaces 
the contact potential g8 3 ('x.)/2 by the pseudo potential 
<7<5 3 (x) (d/dx) x/2 and one can avoid the ultraviolet di- 
vergences by evaluating the partial derivative at the right 
stage of the calculation. We shall not discuss this method 
any further. 

Using Eqs. (53} , (|8*T|) and (|A.7|I . we obtain the num- 
ber density in the one-loop approximation 



M 



d 3 k 



2 ./ (2tt)3 
" M d 3 k 



k^k 2 + 2/z - k 2 - , 



2fc 2 



(2tt) 3 



2 1 



2fc 2 



AiJf 



(84) 



The first integral is now convergent in the limit M — > 
oo and the divergences have been isolated in the second 
integral. This first term goes like M 5 and is independent 
of pL and g. It can therefore be removed by a vacuum 
energy counterterm A\£. The form of AiT can be found 
by substituting g —y g + Aig and /z — > tz + Ai/z in J-o(fi) 
and expanding to first order in Aig and Ai/z. Including 
the vacuum counterterm, we obtain 



«o+i (p) 
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a 



3tt 2 



(90) 



Inverting Eq. (|9()Jl to obtain the chemical potential as a 
function of the number density, one finds 



Mo+i(") = 9^ + -gh.i(2gn) 



Swan 



32 /net 3 
~3 



(91) 



A^(/i) = -ZAw+g-rAig + AiE 
9 V 



(85) 



The counterterms needed to cancel the quintic, cubic and 
linear divergences can then be found by inspection. One 
obtains 



Aig 
Aim 



1 f M d 3 k 2 


(86) 


2 J (2vf k ' 


1 2 f M d 3 fc 1 


(87) 


2 5 J (2tt) 3 fc 2 ' 


1 f M d 3 k 




2 9 J (2tt) 3 ■ 


(88) 



Note that the counterterm Aig in Eq. I|87|l is precisely 
what is needed to cancel the divergence appearing in the 
one-loop correction to the scattering length, Eq. (|30|l . 
This is in accord with our comment on counterterms in 



Note that here and in the following, we are replacing the 
argument \i by gn in the loop integrals I m .n, since correc- 
tions are of higher order in V na 3 . Using Eqs. I|54|) . I|81|) 
and (|90|l . we obtain the energy density in the one-loop 
approximation: 



£o+i(n) 



-I Q -x{2gn) 



= - 2 9n 



= Allan 



128 

Is" 



7T 



(92) 



The leading quantum correction to the mean-field 
result (EU) was first de r ived by Lee, Huang and 
Yang l)Lee and YaneL Il957t iLee Huang and Yang! 11957]) 
for a hard-sphere potential. Later it has been shown that 
it is universal in the sense that it appli es to all short- 
range potential with scattering length a ijBeliaevl Il958t 
iBrueckner and Sawadalll957tlLierlll963|) . Note that the 
result (|92|) is nonanalytic in the scattering length a. This 



13 



shows that the result is nonperturbative from the point of 
view of "naive" perturbation theory, where one uses free 
particle propagators. It corresponds to the summation of 
an infinite set of one-loop diagrams, with repeated inser- 
tions of the operator ^gv 2 ipf. The structure of these dia- 
grams is the same as the ring diagrams first discussed by 
Gell -Mann and Brueckner for the non relativistic electron 
gas l)Gell-Mann and Brueckner! Il957f) , and are summed 
in the same manner. It is interesting to note that each 
diagram in the series is increasingly divergent in the in- 
frared, but the sum is infrared convergent. 

Using (|37|l and the parametrization i|38|) , the expres- 
sion for the number density becomes 



n = v 2 + -(iPi+^) 



(93) 



The condensate density no is given by the expectation 
value v 2 = | (tp) | 2 At the mean- field level, one neglects the 
fluctuations of an d ^2 and replaces (ip*ip) by | ( V^) 1 2 - 
The total number density is then equal to the condensate 
density. Taking quantum fluctuations into account, this 
is no longer the case. Due to interactions, some of the 
particles are kicked out of the condensate and are not 
in the k = momentum state. The difference n — no 
is called the depletion of the conden sate and is propor- 
tional to the diluteness parameter V na 3 . The one-loop 
diagrams that contribute to the expectation values (ip 2 ) 
and (V'f) are shown in Fig.QJ The solid line denotes the 
diagonal propagator for tpi and the dashed line the diag- 
onal propagator for ip2 ■ The blobs denote an insertion of 
the operator ipf or ip 2 . Taking these one-loop effects into 
account, Eq. l|9"3*)l reduces to 



1 f duj [ d d k 
n 0+1 (n ) = n + -zj -J ^ 



\k)/k 2 



e 2 (k) 



no 



+ -[I hl (2gn Q )+ J_i,_! (2gn )} . (94) 



Using the expressions for the integrals in the appendix, 
we obtain 



n Q+1 (n a ) = n Q 



1+ 3 



now' 



(95) 



The result 195|) for the depletion was first obtained by 
Bogoliubov in 1947 llBoeoliu bovl Il94 7j) . In recent exper- 
iments, Cornish et al. ijCorrilsheT'aU boOOj) were able to 
vary the s-wave scattering length a for 85 Rb atoms over a 
large range by applying a strong external magnetic field 
and exploiting the existen ce of a Feshbach resonance at 
B ~ 155 G. Values for V na 3 up to approximately 0.1 
were obtained and should be sufficiently large to see de- 
viations from the mean field in experiments. In order 
to observe this quantum phenomenon, it is essential that 
experiments are carried out at sufficiently low tempera- 
ture so that the thermal depletion of the condensate is 
negligible. 




FIG. 4 One-loop diagrams contributing to expectation values 
(tf) and (iPl). 



D. Collective Excitations 

The Bogoliubov spectrum is given by Eq. I|tj7|) and 
was derived from the microscopic theory represented by 
the action (j26(l . It is linear for small momentum p with 
slope \f2jl. The spectrum uj(p) of collective excitations 
is given by the poles of the propagator. The poles are 
the solutions to the equation 



det [^ - 1 (^, P )-n( w ,p)] 







(96) 



where Do(uJ,p) is the real-time version of the free prop- 
agator 10 and H(uj,p) is the 2x2 self-energy matrix. 
The dispersion relation uj(p) is generally complex and 
can therefore be written as 



uj(p) = Rew(p) — ijip) 



(97) 



The real part Re tu(p) gives the energies of the excitations, 
while the imaginary part "f(p) represents the damping of 
the excitations. The Bogoliubov spectrum is purely real 
and in this approximation, the excitations therefore have 
an infinite lifetime. 

In the following, we calculate the leading quantum 
correction to the real part of the spectrum in the 
long-waveleng th limit and thus reproduce Beliaev's 
classic result (Bcliaev, 1958|) . The one-loop Feynman 
diagrams that contribute to the self-energies are shown 
in Figs. ISlEf. The solid line denotes the diagonal 
propagator for ipi and the dashed line the diagonal 
propagator for ip2- The off-diagonal propagators for ipi 
and ip2 are represented by lines that are half solid and 
half dashed and vice versa. One-loop contributions to 



self-energies are down by a factor of 



compared 



to the mean-field terms in the inverse propagator. It is 
therefore consistent to evaluate the self-energies using 
the mean-field dispersion relation e(p) since corrections 
would be suppressed by at least a factor of na 3 . 



Note that the sum of all one-particle reducible diagrams that con- 
tribute to the self-energies vanishes. This sum is proportional to 
the derivative of the effective potential and is zero when evalu- 
ated at the minimum of the effective potential. 
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FIG. 5 One- loop diagrams contributing to the self-energy 

nii(w,p). 



FIG. 6 One- loop diagrams contributing to the self-energy 

U 2 2(u,p). 



FIG. 7 One- loop diagrams contributing to the self-energy 
IL 12 (uj,p). 



In the following, we calculate the off-diagonal self- 
energy Ui 2 (co,p) explicitly. The expression is 



n 12 (e(p),p) 



2 2 f duj f d d k 
3 [lj + e(p)] k 2 



1 [(lj + e(p)) 2 - e 2 (|p + k|) + ie}[u 2 - e 2 (fc) + ie] 

[Lj + e(p)}e 2 (k)/k 2 1 

[(uj + e(p)) 2 - e 2 (|p + k|) + is] [lj 2 - e 2 (k) + ie] | 



(98) 



After integrating over lj, the imaginary part of 
IIi2(e(p),p) becomes 



ImU 12 (e(p),p) 



2 2 
"3 V 



d d k 



+ 



(2w) d 

e(p) + e(k)][3k 2 /e(k)-e(k)/k 2 ] 
[e 2 (\p + k[)-e 2 (p)-e 2 (k)} 

3k 2 -e 2 (k)/k 2 \ 



[(e(k)+e(p)) 2 -e 2 (\p + k\ 



(99) 



Finally, we expand Eq. (|99|l in powers of the external 
momentum p. After simplifying using Eqs. |A~7|) (TOL 



the self-energy reduces to 



lmU 12 (e(p),p) = i g (2. gt ; 2 ) 3 / 2 [3/i, 3 - /_ M ]; 



+0 (p 3 /g 3/2 v 3 



(100) 



Notice in particular that II 12 (0,0) vanishes. This prop- 
erty holds to all orders in perturbation theory and follows 
from time-reversal invariance. Note also that n 21 (u>,p) = 
—U 12 (lj,p), which also holds to all orders in loop expan- 
sion. The real part of the self-energies Hu(e(p),p) and 
Il22(e(p),p) are expanded about zero external momentum 
in the same way. Including the mean-field self-energies, 
one finds: 



ReU n (e(p),p) = 3gv 2 



1 

A 9 



3/ M + - gv 2 (9J a ,a 



ReH 22 (e(p),p) 



-6/ 0> i + /_ 2 ,-i) 
o 1 



gv 



12 



+ 0(p 2 /gv 2 ) (,101) 
9/i.i + 3/-i,_i 

i L 

-pV 2 (6/3,5 - H/1,3 + 3/_i,i) 
+0(p i /(g 2 v i ) . (102) 



The expressions for the self-energies Hn(e(p),p), 
Il22(e(p))P) an d IIi2(e(p),p) are infrared divergent. The 
integrals /_2,-i(2<?no) and /_i,i(2grto) both have a loga- 
rithmic divergence as the loop momentum k goes to zero. 
These divergences show up as poles in e in Eqs. (|A~T6|) - 
(|A.17|) . However, it is important to point out that they 
cancel in the final results for physical quantities, as we 
shall see below. 



The real part of equation (|96|l can now be written 



as 



[ W -Imn 12 ( £ (p),p)] 2 - [p 2 -fi + ReU n (e(p),p)] 



x P 



Ren 22 (e(p),p)] . (103) 



The next step is to eliminate the chemical potential from 
Eq. (|103fl by minimizing the one-loop thermodynamic po- 
tential Qo + i(/j,,v), This is found by differentiating the 
sum of Eqs. (|73(l and i|80[l with respect to the condensate 
v and setting it to zero. This yields 







-/i + gv 2 + -g 



d d p 3(p 2 -Y) + (p 2 - X) 
(2ir) d ^p^X)(p 2 -Y) 



(104) 



It is consistent to the order in quantum corrections at 
which we are calculating, to evaluate the one-loop con- 
tribution to p04(l at the classical minimum v = vq. 
Eq. I|104|l then reduces to 



= -/i + n 22 (o,o) . 



(105) 
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Eq. (|105[l shows that the Goldstone theorem is satisfied 
at the stationary point of the thermodynamic potential. 
Using 1)94(1 . we see that the chemical potential obtained 
from l|91|) agrees with l|105f) . Solving for the chemical 
potential and substituting the result as well as the self- 
energies given by Eqs. (|100|l - l(l()2|l into Eq. 1(103(1 . we 
obtain 



u + \g{2gv 2 f/ 2 (3/i, 3 - /_i,i)p 



P z + 2gv A ( 1 - ig(9/ 2 , 3 + I- 2 ,-i - 6I ,i] 



l-^V (6/3,5 - ll/i,3 + 3J_i,i; 



(106) 



We can now solve for ui in Eq. i|106|) . In the long- 
wavelength limit, p <C y /, 2gv 2 , we obtain 



Rew(p) = p\/2gv 2 



2 2 

1 - (6/3,5 + 7/i,3 " 3/_i,i) 



1- -:3(9/2.3- 6/0,1 +/-2.-1; 
16 



7\/2gn 



28 n a 3 



3 V 7T 



(107) 



The infrared divergent terms that appear in the expres- 
sions for the self-energies cancel algebraically after hav- 
ing used the relations (|A.7|> — (|A.9I) . The ultraviolet di- 
vergences associated with the integrals in Eq. (|107J) are 
again power divergences. Thus one immediately obtains 
a finite result and this is another example of the conve- 
nience of employing dimensional regulari zation. The re - 
sult (|1U7|I was first obtained by Beliaev (Beliaev, 1958). 
One can easily check that the slope of u> in Eq. (|107f) is the 
same as the macroscopic speed of sound, c, that one ob- 
tains from differ entiating the pressure with respect to the 
number density (Fetter and Waleckalll97l[) . This equal- 
ity has been proven to all orders in perturbation theor y 
by Gavoret and Nozieres llGavoret and No zierest f l964'l. 
and b y Hohenberg and Martin l|Hohenberg and Martini 
119651) . 

The effective field theory approach presented in this 
section is based on a cartesian parametrization of the 
quantum field ip in Eq. 1)38(1 . A similar effective field 
theory approach was develo ped by Popov llPopovl Il987j) 
and later extended by Liu (|LiuL Il997l Il998|) . Instead of 
using a cartesian parametrization of the field ip, it was 
parametrized using the density n and the phase x : 



V>(x,t) = v/^CM) e 4x(x,t) • 



(108) 



The density field is shifted in analogy with the shift in 
Eq. (23): 



n(x, t) = v 2 + er(x, t) 



(109) 



We recall that the condition (|56|l simplifies calculations, 
since it makes all one-particle reducible diagrams vanish. 
This condition is replaced by 



M = 0. 



(110) 



Using the parametrization l|108|) together with (|109fl . the 
action l(26() takes the form 



S\v] 



dt / d d x 



1 t 1 , ■ ■ ; 1 2 
-J<T+7> (X°- XV) - T.9V 
v 2 2 



+ 4>^ 



(111) 



where the dots indicate an infinite series of higher or- 
der operators. The classical action S[v] and the source 
J are given by Eqs. (|4*0|) and (03J), respectively. After a 
rescaling of the fields a and x, the free propagator cor- 
responding to the action ((111(1 is identical to Eq. 1(66(1 . 
However, the interaction vertices are different. In partic- 
ular, the vertex corresponding to the operator cr(Vx) 2 
is momentum dependent. One feature of the pertur- 
bative expansion that follows from the action ((111(1 is 
the ab sence of infrared divergences in indiv idual dia- 
grams l|Braaten and Nietd . ll99alPonovill983l) . The mo- 
mentum dependence of the trilinear interaction er(Vx) 2 
compensates the singular behavior of the propagator at 
low momenta. This difference should not be viewed as 
being fundamental since the infrared divergences always 
cancel in physical quantities. We have already seen one 
example of this, when we considered the quantum correc- 
tion to the Bogoliubov dispersion relation. On the other 
hand, individual diagrams are more severely ultraviolet 
divergent, but these cancel when the diagrams are added. 
The above features merely represent a different way of 
organizing the perturbative calculations. The equality 
of calculations of physical quantities order by order in 
perturbation theory simply reflects the reparametrization 
invariance of the functional integral. 

As mentioned before, the imaginary part of the dis- 
persion relation represents the damping of the collective 
excitations. In the case of a dil ute Bose gas, this term 
was first calculated by Beliaev ( Be liaevt Il958|) . In the 
long- wavelength limit, his calculations showed that the 
damping rate is proportional to p 5 : 



lip) 



3p 5 



3207mo 



(112) 



The imaginary part 7(p) is connected with one phonon 
decaying into two phonons with lower energy. It is of- 
ten referred to as Beliaev damping. The action 11111) 
has later bee n used by several authors ijLiul 119971 1998: 
iPopovt Il983(l to rederive the result |{TTg|| . The calcula- 
tions represent a significant simplification compared to 
the original derivation. 
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E. Nonuniversal Effects 

In the previous subsection, it was shown that the 
dominant effects of the interaction between the atoms 
could be subsumed in a single coupling constant called 
the s-wave scattering length. Thus all interatomic po- 
tentials with the same s-wave scattering length will have 
the same properties to leading order in the low-density 
expansion, and this property is called universality. How- 
ever, at higher orders in the low-density expansion, phys- 
ical quantities will depend on the details of the inter- 
atomic potential such as the effective range r s . These are 
called nonuniversal effects. A detailed analysis of nonuni- 
versal effects can be fo und in the paper by Braaten, 
Hammer and Hermans l|Braaten et all l200l[) . We dis- 
cuss these next. 

Including the operator [V( , 0*'0)] 2 m Eq. J2HJ), we 
can again calculate exactly the scattering amplitude for 
s-wave scattering. Summing t he contributions from the 
diagrams in Fig. we obtain ijBraaten et all l200lfl 



T(q) 



1 ■ g 

2g + Ahq 2 % I6n 



(113) 



The coupling constant h is th en related to the effe ctive 
range r s of the true potential ijBraaten et all 1200 l|) and 
it is determined by matching Eqs. (|34|1 and (|113|l through 
third order in q. This yields 



h = 2na 2 



(114) 



After performing the shift (|37|l . the operator [V(ijj*ip)) 2 
also contributes to the free part of the action. Includ- 
ing the effects of this operator, one obtains a modified 
propagator and a modified dispersion relation: 



IU1 



e 2 (p) +ie V iuj e2 (p)/P 
where the new dispersion relation is 

e(p) = p^(l + 2hv 2 )p 2 + 2gv 2 . 



(115) 



(116) 



The spectrum (|116[) has the Bogoliubov form with mod- 
ified coefficients. It is therefore straightforward to re- 
calculate the ground state energy density and the de- 
pletion of the condensate by rescaling the momentum 
p — * p-\/l + 2hv 2 . For instance, the expression for the 
number density is in analogy with equation (|94|l given by 



, \ 1 
?i +i (no J = n a + - 



d d p p 2 e 2 (p) 



{2ir) d e(p) p 2 



n + X - [1 + 2hv 2 ] ^ I 1A 



[l + 2hv 2 



I-i, 



(117) 



The limit d — > 3 is regular and we obtain 
«o+i (no) = n a 



1 y/8n g 3 , 2 x 
1+ 24^ (l + 2^) 2 ( 1 - 2 ^) 



(118) 



Using Eqs. 1|36|) and (f 1 1 4|) . and expanding to first order 
in the effective range r s , we find 



"o+i(™o) = n 



n ^ 

3 V 7T 



32ir 2 r s ( n§a 



3 \ 3/2' 



(119) 



Similarly, we can calculate the ground state energy den- 
sity. Exp anding to first order in the effective range r s , 
one finds ijBraaten et all l200lj) : 



£o+i(n) = Airan 2 



1024n 2 r s 
15o 



128 / 


na 3 




7T 


^na 3 


^3/2" 







(120) 



Effective field theory can be used to determine at which 
order in the low-density expansion a given operator starts 
to contribute to various physical quantities. As an exam- 
ple, we consider the energy density £. Each power of 
contributes a factor of n. Each power of V contributes 
a factor of \fna~. Each loop order in the quantum loop 
expansion contribitutes a factor of V na 3 . The derivative 
interaction [V(ip*ip)] 2 does not contribute to the energy 
density at the mean-field level for a homogeneous Bose 
gas since Vv obviously vanishes in this case. It first con- 
tributes^ at the one-loop level and this gives one factor 



There are two powers of ip*ifj and two powers 



of V, which give two factors of n and ^/na, respectively. 
This yields a contribution proportional to r s n 2 (na 3 ) 3 / 2 
in accordance with the explicit calculation (|120H • 

We next comment on the nonuniversal effects that 
arise in higher-order calculations of the energy density. 
For simplicity, we ignore all other coupling constants 
than g. The leading term is the mean-field contribu- 
tion given in (|76|l . A one- loop calculation gives rise to a 
universal correction proportional to V na 3 shown in 192fl . 
At the two-loop level, one encounters two terms. The 
first is a universal logarithmic corrections proportional to 
na 3 lo g na 3 . This term was first calculated by Wu (Wu, 
11959ft . The second term is a nonuniversal term propor- 
tional to n a 3 . This term was first ca lculated by Braaten 
and Nieto ( Braaten and Nieta .il 997). It comes about as 
follows. At the two-loop level, there is a logarithmic ul- 
traviolet divergence that cannot be cancelled by a local 
two-body counterterm of the form Sg(ip*ip) 2 . The only 
way to cancel it, is to add to the Lagrangian a local coun- 
terterm of the form 5g^,(ip*ip) 3 and absorb the divergence 
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in the coefficient of this operator. One can see the neces- 
sity of such an operator by considering 3 — ► 3 scattering. 
At the two-loop level, or fourth order in a, there are 
Feynman diagrams that depend logarithmically on the 
ultraviolet cutoff. They give an additional momentum- 
independent contribution to the 3 — ► 3 scattering ampli- 
tude. In order to reproduce the low-energy scattering of 
three atoms, the local momentum-independent operator 
Sg^faip) 3 must be included in the effective Lagrangian. 
The counterterm of this operator that removes the loga- 
rithimic divergences from the 3 — > 3 scattering amplitude 
is then exactly the same counterterm needed to remove 
the logarithmic divergence in the energy density. 

The coefficient g% generally depends on the proper- 
ties of the two-body and three-body potentials. The op- 
erator (fafa 3 m Eq. <|26|) takes into account not only the 
contribution from 3 — > 3 scattering from a possible three- 
body potential, but also the contribution from the succes- 
sive 2 — > 2 scattering via the potential Vo (x) . One way to 
determine the coefficient g$ would be to solve the 3 — > 3 
scattering problem for the potential Vo(x). Alterna- 
tively, one can determine it from calculating the ground 
state energy density of bosons interacting through Vo (x) . 
Such a strategy was recently used by Braaten, Ham- 
mer and Hermans ijBraaten et flfl. 120011) using the Monte 
Carlo calculations of the condensate fraction and energy 
density for four differen t model potentials by Giorgini, 
Boronat, and Casulleras IjGiorgini et a/l ll999). The four 
potential were a hard-sphere potential with radius a, two 
soft-sphere potentials with height Vo and radii R = 5a 
and R = 10a, and a hard-sphere square-well potential 
with depth Vo with inner and outer radii of R = a/50 
and R — a/10, respectively. These four potentials all 
have the same s-wave scattering length a, but different 
effective range r s . By calculating the energy density for 
the homogeneous Bose in the low-density expansion and 
matching it onto the Monte Carlo results, Braaten, Ham- 
mer and Hermans were able to estimate the coefficient 
gs. Due to the large statistical errors, they could not 
find any deviation from universality in the three-body 
contact parameter. In order to determine the coefficient 
more accurately, one needs data with higher statistics at 
various densities. 



A. Hartree-Fock Bogoliubov Approximation 

The self-consistent Hartree-Fock-Bogoliubov ap- 
proximation and its relation to the Bogoliubov approx- 
imation and the Popov approximation have been dis- 
cussed in detai l jGriffinl Il996t iHutchinson et alj l2000t 
IShi and Griffiri Il998|) . The starting point is the ac- 
tion (|26|) in imaginary time: 



S[fa,ip] = dr d d x\ ip 



9 rv2 



(121) 



The next step is to treat the interaction term using a 
self-consistent quadratic approximation. After perform- 
ing the shift l|37l) . there are terms that are cubic and 
quartic in the quantum field fa These terms are approx- 
imated as follows 

fa^ w 2<-0*-0>-0 + (-0-0)^* , (122) 
fatpfaip w 4:{tp*lp)tp*lp + (fafa)ipip + (tpfafafa . 

(123) 

Terms involving the expectation value of a single field 
have been omitted since (fa — (fa) — 0. Moreover, 
terms involving the expectation values of three or four 
fields are omitted. The quantities (fa fa and (ipfa are 
often referred to as the normal and anomalous average, 
respectively. 

We next write the quantum field tp — (fa ~ s rifa)/V^ 
and insert Eqs. I|122(l - (|123|) into the action !|121|) . One 
then obtains an approximate action which is quadratic 
in the fluctuating fields: 

S[v,lpi,1p 2 ] = S[v] + Sfree[v,1pl, Ips] + Si nt [v , Ipi , fa] , 

(124) 

where in analogy with the zero-temperature case, we have 
defined the classical, free, and interacting parts of the 
action by 



IV. WEAKLY INTERACTING BOSE GAS AT FINITE 
TEMPERATURE 



In this section, we discuss the weakly interacting 
Bose gas at finite temperature. We first review the 
the Hartree-Fock-Bogoliubov (HFB) approximation, the 
Bogoliubov approximation, and the Popov approxima- 
tion. We then discuss Wilson's renormalization group 
approach applied to this problem. Finally, we discuss 
improved variational approaches to the finite tempera- 
ture Bose gas. 



S[v] = J* ' dr J d d 



S hee [v, ipifa] 



dr / d d x 



1 A 

2 y 



2 



dr / d d x< 



fiV 

fa (-V 2 + X) <tp 



(125) 



+ 2^2 (-V 2 + Y) fa 



(126) 



H + g(v 2 + 2(fafa 
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\f2vfa 



1 



(fa fa) - y/2wk } , (127) 



where 



W = 
X = - f i + g 

Y = -fi + g 



1 , r 



v 2 + 2(fafa - -$fa - -(fa fa) 



(128) 
(129) 

(130) 



Demanding that the terms linear in fa and fa in 
Eq. i|127|) vanish, immediately yields 



o = (rr) - (w) , 



= —fi + g 



(131) 
(132) 



The equation of motion for the quantum field ip follows 
from (|126|) and reads 



f r = -[v 2 + ^ + ff 



(v 2 + (fa fa) ^ 



+ ^ + {i>ip))ip 



(133) 



The propagator that corresponds to the free part of the 
action is 



D(uJ n ,p) 



1 



p 2 - 2g(^>fa uj r , 



(134) 



where we have used Eq. (|131fl to simplify. The dispersion 
relation is given by 



(p) = p 2 -2«?(^) [p 2 + 2 ff « 



(135) 



It can be shown tha t the anomalous average (ipip) is neg- 
ative jGriffml . 11996). so that the dispersion relation 1|135|) 
makes sense. 

The Hartree-Fock-Bogoliubov approximation is 
given by Eqs. I|132|) (|135H . One of the attractive fea- 
tures of the HFB approximation is that it is a conserving 
approximation that is guaranteed to respect the conser- 
vation laws that follow from the underlying symmetries 
of the field theory. One of the problems with the HFB ap- 
proximation is that there is a gap in the spectrum (I135fl . 
which is a consequence of the spontaneously broken sym- 
metry. Thus the HFB approxima tion violates the Gold- 
stone theorem. Griffin and Shi ijShi and Oriffirl 11998ft 



have argued that the problem of a gap in the HFB ap- 
proximation is due to overcounting diagrams that con- 
tribute to the anomalous self-energy and are second or- 
der in the interaction. Thus the HFB approximation is 
inconsistent to that order in the interaction. 

Another problem with the HFB approximation is 
that is it computationally difficult to apply. One starts 
with some initial guess for the condensate density, the 
chemical potential as well as the normal and anomalous 
averages satisfying Eq. I|132|) . One then solves for the 
infinitely many normal modes in Eq. I|133|) . These nor- 
mal modes are then used to calculate the normal and 
anomalous averages. The procedure is iterated to self- 
consistency d . In this way, one obtains information about 
the different normal modes. However, since the HFB ap- 
proximation has a gap, some of this information must be 
qualitatively incorrect. 

The HFB approximation can also b e viewed as a 
variational one l)Blaizot and Ripkal IT986). although the 
equivalence seems to have gone unnoticed in some of the 
literature. The idea is to make a Gaussian ansatz for 
the ground state wave functional and the excitations. 
The variational parameters are the normal and anoma- 
lous self-energies En(w„,p) and Ei2(w„,p). They are 
normally taken to be independent of frequency and mo- 
mentum. The corresponding terms Hn(uj n ,p)faij} and 
Si2(w n) p)V'' ! / ; are quadratic in fields and are simply added 
to and subtracted from the Lagrangian. Expressing the 
normal and anomalous self-energies in terms of tin and 
H22 , the action can be split into a free part and an inter- 
acting part according to 



5 , f ree [v,t/'iV'2] = / dr I d d x 







-i (^fafa - ^1^2) 



~fa (-V 2 - X) fa + ~fa (-V 2 - Y) fa 



, (136) 



Sud,[v,fa, fa] = I dr I d d x 



- (3gv 2 - nn) fa 2 



1 (gv 2 - n 22 ) fa - ±=Zfa (fa + fa) 



1 



where we have defined 



X = /! 
Y = n 

Z = -gv 



Tin 

n 22 

2 



(137) 



(138) 
(139) 
(140) 



d In practice one faces the problem of ultraviolet divergences in 
the anomalous average. Various rcnormalization procedures ha ve 
been discussed by Hutchinson et al. (Hutchinson et al, 2 0001) . 
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The propagator that corresponds to the free part of the 
action is 



D(u n ,p) 



1 



e 2 (p) 



P 



Y uj n 

n P 2 -X 



where the dispersion relation is given by 
e(p) = y/(p*-X)(j?-Y). 



(141) 



The next step is to calculate the thermodynamic poten- 
tial in some approximate way including the interaction 
term 1(137(1 which consists of three and four-point ver- 
tices together with the subtracted self-energies. This 
is done by calculating the thermodynamic potential ap- 
proximately according to 



il = —p,v' 



1 



-gv 



TrlogZT 



(Si 



hit/ 



(142) 






FIG. 8 Vacuum diagrams included in the Gaussian approxi- 
mation. 



dTh~2 



(147) 



where (A) is the thermal average of the operator A. The 
thermodynamic potential then becomes 
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-I 9 
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i^ + iTrlogi}- 1 
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2 (p) 



p 2 -Y 
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P u$ 



t 2 {p) 



:(3gv 2 
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=2(p) 

X 



'■(p) 



(143) 



The Feynman diagrams that correspond to the interac- 
tion term in Eq. i|143|) are shown in Fig. [5] The blob 
denotes an insertion of either 3gv 2 — Iln or gv 2 — H22- 

The condensate density v is determined by the sta- 
tionarity condition 



dv 



(144) 



Using the expression l|143|) for the thermodynamic po- 
tential, Eq. 11441) becomes 







-/! + gv 2 + ^gjj, 
~2 9 4p [J n + e 2 (p) 



p 2 -Y 
>1 + £ 2 (P) 



(145) 



The self-energies Iln and II22 are determined variation- 
ally by demanding that they minimize the free energy: 



0, 



(146) 



The equations ((146(1 and l|147|) arc often referred to as 
gap equations. The solutions are 



nn 



n 22 = 



3gv 2 
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p 2 -Y 



e 2 (p) 



p 2 -X 



^7 



•e»(p) 
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i4 P 



P L^ + e 2 b) 
p 2 -X 



^7 



e»(p) 



(148) 



(149) 



The Feynman diagrams that correspond to the gap equa- 
tions (|148fl and 1149(1 are shown in Fig. Insert- 
ing l|148[) and (|149|) into l|141|) and using (|145|l to elimi- 
nate the chemical potential, we obtain the Hartree-Fock- 
Bogoliubov spectrum 1(135(1 . 



FIG. 9 Feynman diagrams included in the self-energies II11 
and II22 in the Gaussian approximation. 

The Gaussian approximation can be made the start- 
ing point for a systematic expansion procedu re. This ex 
pansi on was first formulated by Okopinska llOkopinsk 
1987). Stancu and Stevenson ijStancu and Stevenso: 
1990) calculated the leading corrections to the Gaussian 
approximation for a relativistic 4 -theory in four dimen- 
sions. The Feynman diagrams that are included in this 
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next-to-leading order calculations are the three-loop di- 
agrams, the two-loop setting sun diagrams, the two-loop 
double bubbles with a single self-energy insertion, as well 
as the one-loop diagrams with two insertions of a self- 
energy. 

We next discuss the Bogoliubov and Popov approx- 
imations, which are obtained by making certain approx- 
imations in the full HFB approximation. 



The other thermodynamic functions follow from the free 
energy 11521) . In the low-temperature limit, the number 
density becomes 



9 



3tt 2 



7T 2 T 4 



30(2 M ) 5 / 2 



(154) 



Inverting (|154fl . we obtain the chemical potential in the 
limit T <C 2gn: 



1. Bogoliubov approximation 

The Bogoliubov appoximation amounts to ne- 
glecting both the normal and anomalous averages in 
Eqs. (|H2V((H5|> . The chemical potential ((HISjl then re- 
duces to 



fi = gv . 

The spectrum (|135|) is given by 



-(p) = pVp 2 + 2 M , 



(150) 



(151) 



which is gapless. The Bogoliubov approximation has a 
natural interpretation in terms of Feynman diagrams. 
It is obtained by neglecting the two-loop diagrams con- 
tributing to the thermodynamic potential l|143|) . The 
self-energies that follow from the gap equations are then 
given by their mean-field values, so we have Un(u> n ,p) — 



3gv 2 and n 22 (w„,p) 



The stationarity condi- 



tion (|144f> simply reduces to i|15U|) The free energy is 
obtained by substituting the mean-field values for the 
self-energies into the thermodynamic potential evaluated 
at the minimum. This yields 



T = 



' Jp 



(152) 



In the Bogoliubov approximation, one assumes that most 
particles are in the zero-momentum state. Clearly, this 
approximation is only valid at very low temperatures 
where one can ignore the thermal depletion of the con- 
densate. 

One can expand the free energy (|152|) about zero 
temperature. At sufficiently low temperatures, the ther- 
modynamics is dominated by the phonon part of the 
spectrum. Using Eqs. HA.10|) and l|A.31fl . the free energy 
reduces in the limit T <C 2n to 



T 



4yW 

15tt 2 



7T 2 T 4 



90(2^t) 3 / 2 



(153) 



M = 9 n 



1 , yZng 
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TT 2 T 4 
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(155) 



Similarly, for temperatures T <C 2gn, the equilibrium 
energy density is 



£ = Airan 2 



128 na 3 



15 V 7T 



tt 2 T 4 



30(2 5 n) 3 /2 



(156) 



Using Eq. (|93|l . one can calculate the total number den- 
sity as function of the condensate density and tempera- 
ture: 



n = n Q + 



P 2 + e 2 (p) 
2% luI + £ 2 {p) 



no + - + 

1 [°° P 2 (P 2 + 6 2 ( P )) 
+ 4^X dp ifo) n{t{p)) 



In the limit T -C 2gno, one finds 

T 2 



(157) 



n = no 



24^2^ 



(158) 



Eqs. (11531) - (11581) w ere fi rst obtained by Lee and 
Yang ijLee and Yand . ll95o|) . The second term inside the 
brackets is the quantum depletion of the condensate that 
was calculated in Sec. IIII.Cl The last term is the thermal 
depletion. 



2. Popov approximation 

In the Popov approximation, one neglects the 
anomalous average in Eqs. I|132|) (|135l) . The chemical 
potential (|132fl reduces to 



M = 9 



v 2 + 2(^>) 



Similarly, the spectrum l|135|l becomes 



~{p) = pVp 2 + 



(159) 



(160) 
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and is gapless. Note that the spectrum formally has the 
same form as the Bogoliubov spectrum, but now the con- 
densate density depends on the temperature. The num- 
ber density n satisfies 

n = n + . (161) 

Eqs. H159|) and l|161|) constitute the equation of state for 
the weakly interacting Bose gas in the Popov approxi- 
mation. For a given temperature T and total number 
density n, they must be solved simultaneously for the 
condensate density no and the chemical potential. Note 
in particular that the condensate density has a strong 
temperature dependence. For small values of the gas pa- 
rameter, the condensate density as a function of temper- 
ature deviates typically only by a few percent from the 
result 118fl for the ideal gas. 

At T = 0, the expectation value in Eq. (|159fl is sup- 
pressed by a factor of V na 3 compared to the mean-field 
term. This implies that the Popov approximation gives 
the same results as the Bogoliubov approximation for 
all thermodynamic q uant ities up to corrections of order 
the gas parameter V na 3 . Since the Popov approxima- 
tion does not include all corrections of order V na 3 , it is 
no more accurate than the Bogoliubov approximation at 
zero temperature. 

The Popov approximation ought to give a good de- 
scription at higher temperatures than the Bogoliubov ap- 
proximation since it takes into account the temperature 
dependence of the condensate. However, the Popov ap- 
proximation breaks down in a narrow temperature region 
around T c . This can be seen from the fact that the Popov 
approximation predicts a first-order phase transition for 
the weakly interacting Bose gas, while universality argu- 
ments based on the 0(2)-symmetry tell us that the phase 
transition is of second order c . 

The Popov approximation can also be interpreted 
in terms of Feynman diagrams. This is done by ex- 
pressing propagators and interactions in terms of normal 
and anomalous self-energies rather than rin and II22- 
The Feynman diagrams that contribute to the thermody- 
namic potential in the HFB approximation are still those 
shown in Fig. [8] but the sum-integrals are now in terms 
of normal and anomalous propagators and the symmetry 
factors are different. The Popov approximation for the 
thermodynamic potential is then defined by keeping the 
diagrams that involve normal propagators and neglecting 
those that involve anomalous ones. The corresponding 



gap equation that follow from the variational principle, 
gives the normal self-energy as the sum of the mean- 
field contribution 2gv 2 and the one-loop tadpole diagram 
that involves the normal propagator. The anomalous self- 
energy is given by the mean-field contri bution qv 2 alone . 
This definition has been given before l|Shi and G riffin. 
1998). Finally, the free energy is obtained by substitut- 
ing the expressions for the self-energies into the thermo- 
dynamic potential evaluated at the minimum. The free 
energy in the Popov approximation has the Bogoliubov 
form p52|l . 



3. Many-body T-matrix and modified Popov approximation 

In Sec. IIII.AI we calculated the exact scattering am- 
plitude for 2 — > 2 scattering in the vacuum. In the limit 
where the external momentum goes to zero, the two-body 
scattering matrix goes to a constant. In the next sec- 
tion, wc will show using renormalization group methods 
that the effective coupling constant for a weakly inter- 
acting Bose gas is temperature dependent, and in par- 
ticular that it vanishes at the critical temperature. This 
behavior is expected at a second-order phase transition 
where the correlation length goes to infinity v(the ef- 
fective chemical potential or the effective mass goes to 
zero) due to the fact that ^-theory is a trivial theory. 
It is therefore important to improve on the Popov the- 
ory by using an effective temperature-dependent coupling 
constant. This can be done by using the many-body 
T-matrix. The many-body T-matrix takes medium ef- 
fects into account by summing repeated two-body scat- 
tering processes of quasi-particles in the gas rather than 
in the vacuum. This is done by calculating the dia- 
grams in Fig. |21 at finite temperature using the prop- 
agator The many-body T-matrix T MB (k, k, K; z) 
depends on the relative momenta of the two atoms k and 
k before and after the collision, and the total center-of- 
mass momentum K and center-of-mass energy z. In the 
following, we neglect this energy and momentum depen- 
dence which is a good approx imation l)Biilsma and Stool 
h99fiaHShi a,nd GriffiniriM. 

At zero external momentum p, the one-loop diagram 
in Fig. is 

= 9h,i +g\l dp^n(e(p)). (162) 



In Sec. IIV.CI we show that the long-distance properties of the 
dilute Bose gas is given by a classical field theory in three di- 
mensions with an 0(2) symmetry. In particular, the infrared 
Wilson-Fisher fixed point is that of a three-dimensional 0(2) 
model. This model is known to have a second order phase tran- 
sition. 



By summing the geometric series co rrresponding to the 
diagrams shown in Fig. |S1 one finds ijBiilsma and Stool 
h996aHShi and GriffinllTO 



T (0) 



2 5 2 



2.9-^(0) 



(163) 
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We next consider the many-body T-matrix at low tem- 
perature. For T <C 2(7no, we obtain 



T MB (0) = g i 1 + 



4tt 2 



6 Vn g 



(164) 



We note in particular that at T = 0, the many - bod y 



T-matrix reduces to g up to corrections of order yna 3 . 

The many-body T-matrix has been used to obtain an 
impr oved approximation from the Gaussian approxima- 
tion l)Biilsma and Stool . ll996a|l . The normal and anoma- 
lous self-energies can be easily calculated from Eqs. 114811 
and JHOJl. We obtain 
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E12 = gv +g 
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The last equation can be easily solved for £12. Using 
(|163[1 . we obtain 

S X2 = T mb (0)t; 2 . (166) 

Thus the normal self-energy £n(0,0) is given in the 
Hartree-Fock or in the one-loop approximation, while 
the anomalous self-energy £12(0, 0) is given in the many- 
body T-matrix approxima tion. This observation m oti- 
vated Bijlsma and Stoof ijBiilsma and Stooa . Il996a|) to 
redo the calculation using the many-body T-matrix as 
an effective interaction. This gives both the normal and 
anomalous self-energies in the many-body T-matrix ap- 
proximation: 



£ u = 2T MB (0) 
Eia = T MB (0)t> ; 



(167) 
(168) 



These modified self-energies define a modified Gaussian 
approximation, which was used to investigate the ther- 
modynamic properties of the homogeneous Bose gas in 
two and three dimensions. In three dimensions, this 
approach yields a second order phase transition, but 
the critical temperature is the same as that of an ideal 
gas. At the same time, the effective coupling constant, 
which is precisely the many-body T- matrix, vanishes 
at th e temperature given by Eq. JEJ) |£h i and G riffin , 
1998). However, it turns out that the Hugenholz-Pines 
theorem is not always satisfied. At very l ow tem pera- 
ture, it can be shown (Biil sma and Stool Il996al) that 
the value of v 2 that minimizes the thermodynamic po- 
tential does not exactly correspond to the condition 
M = £n(0,0) -£i2(0,0). 



The self-energies (|167|1 and (|168ll can also be ob- 
tained from the HFB equations l|129|) and 11301) by ne- 
glecting the anomalous average and replacing the cou- 
pling constant g by the many-body T-matrix. By mak- 
ing the substituting g — > T MB and neglecting the anoma- 
lous average in the remaining equations that define the 
HFB a pproximation, a gapless approximation was o b- 
tained ((Hutchinson et a/.Ul998tlProukakis et aZJ . 11998) f 
Eqs. p32|) . 1| 133(1 and 1|135|) now become 



= -fi + T 



MB 



v 2 + 2(^*ip} 



^ = _[v 2 + m ]t]; + t mb 

ot 1 J 

+V 2 V] , 
e 2 (p) = p 2 [p 2 + 2T MI V] . 
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(169) 

(170) 
(171) 



The modified mean-field approximation does not repro- 
duce known perturbative results at T — 0. For instance, 
the prediction for the correction to the Bogoliubov spec- 
trum at long wavelengths differs from Eq. i|107|) • Thus at 
T = 0, this approximation is no more accurate than the 
Bogoliubov or Popov approximations. 



B. Other Variational Approaches 

In this subsection, we discuss other variational ap- 
proaches to the weakly interacting Bose gas. The idea 
is to define a thermodynamic potential fl that depends 
on a set of variational parameters a^. The free energy 
and other thermodynamic variables are then given by 
the thermodynamic potential and its derivative evaluated 
at the variational minimum dQ/dai = 0. A variational 
method can be successful only if the essential physics can 
be captured by the variational parameters. 



1. <l>-derivable approach 

The ^-derivable approximation is an approach in 
which the full propagator serves as an infinite set of 
varitional paramet ers. It was first formulated by Lut- 
tinger and W a rd llLuttineer and Ward Il960j) . and by 
Baym ilBavml [1962) for nonrclativistic fcrmions, and 
later generalized to relativist ic field theories by C orn- 
wall, Jackiw, and Tomboulis l|Cornwall et all 119741) . A 
property of the ^-derivable approximation is that it is 



It is not gapless in the sense that the value of the conden- 
sate v that minimizes the effective potential coincides with 
\x = £ii(0, 0) — £12(0, 0), as explained above. 
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conserving, which means that it respects the conserva- 
tion laws that follow from the global symmetries of the 
system. 

We next apply the ^-derivable approach to the di- 
lute Bose gas. We will rederive some of the results 
obtai ned by Lundh and Rammer ijLundh and Rammeil 
2002), but the formulation is somewhat different. 

The ^-derivable thermodynamic potential tt[D] has 
the form 

n[D] = -/iv 2 + ^gv A + ilrlogLr 1 

-^TrlLD + $[£>] , (172) 

where IT is the exact self-energy, D is the exact propaga- 
tor, and the interaction potential $>[D] is the sum of all 
two-particle irreducible (2PI) vacuum diagrams. These 
are diagrams that do not fall apart by cutting two propa- 
gator lines. It is understood that both II and D are 2x2 
matrices. Tr denotes the trace in configuration space. 
The condensate v is found by minimizing the thermody- 
namic potential in the usual way 




FIG. 10 Vacuum diagrams contributing to the one-loop $- 
derivable approximation $i[_D]. 



are extremely difficult to solve. In certain simple cases, 
where the self-energies are momentum independent, one 
can solve these equations. 

Consider the one-loop ^-derivable approximation. 
The diagrams that contribute to ^[-D] are shown in 
Fig. and it reads 

2 Tpw„+£ 2 (p) 2 J P u)l + e 2 {p) 

(177) 

The self-energies are given by 



dil[D] 
dv 



. 



(173) 



Similarly, the variational principle requires that the ther- 
modynamic potential be stationary under variations of 
the full propagator at fixed Dq. This can be expressed 
as 



dn[D] 
dD 



= . 



(174) 



If we denote the free propagator by Do, the Schwinger- 
Dyson equation for the exact propagator D can be writ- 
ten as 



D- 1 = D n 



n 



Using Eq. (|175|l , we can rewrite Eq. (|174|) 



d$[D] 
dD 



(175) 



(176) 



Eq. (|176|) for the self-energy cannot be solved exactly, 
but one can resort to some systematic approximation. 
The n-loop <&-derivable approximation is such an ap- 
proximation. It is defined by keeping all 2PI diagrams 
up to n loops in the thermodynamic potential. Differ- 
entiation with respect to the components of the propa- 
gator is equivalent to cutting the corresponding lines in 
the Feynman diagrams. Thus the gap equations for the 
self-energy in the n-loop ^-derivable approximation con- 
tain all 2PI diagrams up to n-1 loops. It leads to a set 
of integral equations for the self-energies that generally 



n n 

n 22 



3gv 2 , 
3gv 2 , 



(178) 
(179) 



while FIi2 and II21 both vanish. Thus the self-energies are 
those of the Bogoliubov approximation. We next substi- 
tute the self-energies 111780 and H 17911 into the thermody- 
namic potential (|172|) . It turns out that the two terms 
— iTrlLD and ^[.D] cancel each other and the one-loop 
^-derivable approximation for the free energy reduces to 
the Bogoliubov approximation l|152(l . 

The two-loop ^-derivable approximation $2 is very 
complicated. The two-loop diagrams that contribute to 
the thermodynamic potential are shown in Fig. 1111 The 
corresponding equation for the self-energy matrix 1)176(1 
is obtained by cutting the lines in the diagrams. The 
diagrams that are first order in the interaction are mo- 
mentum independent and can be easily calculated. The 
diagrams that are second order in the interaction are 
momentum dependent. The difficult momentum depen- 
dence comes from the self-energy in the propagators. 
This leads to intractable integral equations for the self- 
energies. There have been attempts to simplify the two- 
loop ^-derivable approximation by making the ansatz 
that t he self-energy is a momentum- independent mass 
term ( Am elino-Camelia and Pil Il993|) . However, such 
approximations lead to a mass gap and hence violate the 
Hugenholz-Pines theorem. 

Finally, wc notice that if one neglects the setting 
sun diagrams in Fig. ^] the two-loop ^-derivable ap- 
proximation reduces to the Gaussian approximation that 
was discussed in detail in the previous subsection. 
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FIG. 11 Vacuum diagrams contributing to the two-loop $- 
derivable approximation $2[-D]. 



2. Optimized perturbation theory 

The complexity of the ^-derivable approach calls for 
the use of simpler variational approaches. One introduces 
a finite number of variational parameters. In its simplest 
form, one introduces a single variational parameter m or 
(i. In a relativistic field theory, m is a variational mass 
parameter, while in a nonrelativistic field theory, ji is a 
variation al chemical pot ential. It was first formulated by 
Yukalov The method can be extended 

to include variational cou pling; constan ts and a varia- 
tional kinetic energy term IjChikul 12000). At this point, 
it is important to emphasize that the parameters intro- 
duced in optimized perturbation theory are completely 
arbitrary and that one needs a prescription for them, 
in order to complete a calculation. One prescription is 
the principle of minimal sensitivity (PMS), where one 
requires that the parameters satisfy a stationarity condi- 
tion. For instance, one can demand that the free energy 
be stationary with respect to variations of these param- 
eters. Another criterion is the principle of fastest appar- 
ent convergence (FAC). This condition requires that the 
difference between a physical quantity calculated at two 
different loop orders be as small as possible. 

The starting point is the action l|121|) . We rewrite 
the corresponding Lagrangian by introducing an effective 
chemical potential fj,%, an effective coupling constant A, 
and a counting variable S in the following manner: 



C = V;* 



/ii - 6((i - hi) 



+ Iaoa>) 2 + ^%-A)W>>) 2 + 



(180) 



If we set 5=1, we immediately recover the Lagrangian 
that corresponds to the action (|121|l . Optimized pertur- 
bation theory is defined by the power counting rule that 5 
be of the order A ~ g. We carry out calculations in pow- 



ers of A and at the end of the calculations, we set 6=1. 
It is important to note that the power counting rules 
also should be applied to the counterterms. The ultra- 
violet divergences that appear in optimized perturbation 
theory are removed by the counterterms determined in 
the standard loop expansion for perturbation theory at 
T = 0. 

The introduction of fii and A represents a reorga- 
nization of the perturbative series, which is a selective 
resummation of higher order graphs. But it is important 
to emphasize that there is no overcounting of diagrams. 
Every diagram is counted once with the correct symmetry 
factor. If calculated to all orders, the results for physi- 
cal quantities would be independent of these parameters. 
However, at any finite order in perturbation theory they 
do depend on these parameters. 

Performing the shift l|37|l . the free propagator now 
takes the form 



D(uJ n ,P) 



where 



1 



2 (P) 



p 2 -Y uj n 
-u n p 2 - X 



X = fii — 3Xv 2 , 
Y = \v 2 , 



e(p) = V(P 2 ~X)(p 2 -Y) . 

The mean-field thermodynamic potential is 



-fllV 



-Xv 4 



(181) 



(182) 
(183) 
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(185) 



The one-loop contribution to the thermodynamic poten- 
tial reads 



fii - -5{H-Hi)v 2 + l -8{g-\)v 4 



i^f log [u,l + (p 2 -X)(p 2 -Y)] 



(186) 



The one-loop thermodynamic potential is then given by 
the sum of Eqs. 1)185(1 and 1)186(1 . Setting 6 = 1, we obtain 
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0+1 



-jJLV 
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+i^f log [oj 2 n + (p 2 -X){p 2 -Y)\ 



(187) 



At this point, we would like to discuss the Goldstone 
theorem in connection with optimized perturbation the- 
ory. There is some confusion in the literature whether 
OPT violates Goldstone's theorem. Depending on the 
choice of the parameters fii an d A, the mean-field dis- 
persion relation e(p) in Eq. (|184|l may or may not be 
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gapless when evaluated at the minimum of the effective 
potential. Beyond the mean-field approximation, the dis- 
persion relation is, however, not given by Eq. 1184(1 . but 
rather by 



det [£> - 1 K,p)-n(w„,p)] = 



(188) 



where Dq(uj,p) is the propagator © and H(u) n ,p) is the 
2x2 self-energy matrix. The solution to Eq. I(188[l is 
gapless order by order in optimized perturbation theory. 
We can easily check that at the one-loop level. Differen- 
tiating the effective potential (|187fl with respect to the 
condensate and demanding that v be a stationary point 
yields 



-fj.- 



n 22 (o,o) , 



3{p 2 -Y) + {p 2 - X) 
u>% + e 2 (p) 



(189) 



where 1122(0, 0) is imaginary-time version of (|102fl and 
corresponds to the self-energy diagrams displayed in 
Fig. [fj] Eq. (|189f) ensures that Goldstone's theorem is 
satisfied at one loop. 

We next discuss various choices of the parameters 
Hi and A. Although they in principle are arbitrary, some 
choices are better motivated from a physics point of view 
than others. One requirement one might impose is that 
the results reduce in the limit T — > to the results one ob- 
tains by applying the perturbative framework discussed 
in the previous section. One very simple choice is 



gas at T = 0. Furthermore, it also takes into account 
medium effects by summing repeated two-particle scat- 
tering in the gas. 

Relativistic ^ 4 -theory at finite temperatu r e has 
been s tudied using OPT l|Chiku and HatsudaL Il998t 
IChikul l2000l). The one- loop calcu lation that was car- 
ried out ijChiku and Hatsudalll998j) predicts a fir st-order 
phase transition, while the two-loop calculation (Chiku, 
2000) predicts a second-order phase transition. Thus it 
is very likely that a two-loop calculation using optimized 
perturbation theory is capable of describing correctly this 
important aspect of the phase transition, while incorpo- 
rating the Goldstone theorem. Secondly, OPT provides 
an ideal framework for calculating the energy shifts and 
damping rates of collective excitations of a Bose-Einstein 
condensate. Clearly, however, more work is needed. 

In the paper by Haugset et al. ijHaueset et aZl ll9981. 
the authors construct an improved one-loop thermody- 
namic potential by including the self-energy 1122(0, 0) in 
the propagator. The calculation of the improved one- 
loop effective potential is equivalent to summing up all 
the ring diagrams, except that the two-loop diagram is 
counted twice. One must then subtract by hand the 
term that was overincluded. From the viewpoint of OPT, 
this overcounting comes about because the one-loop di- 
agram with a self-energy insertion was omitted (in ad- 
dition to the setting sun diagrams). Thus a consistent 
power counting according to the rules above is neces- 
sary to avoid problems with overcounting of Feynman 
diagrams. 
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This choice leads to the dispersion relation \67}. Another 
choice is motivated by the principle of minimal sensitiv- 
ity where one demands that the parameters satisfy 



o— =0 ' 

^ = 0. 



(192) 
(193) 



At the one-loop level, the only solution to Eqs. 
(fT35)l is the trivial solution m = A = <|Chikut 12000). 
Only at higher orders are there nontrivial solutions. A 
very promising choice is 
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(194) 
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The many-body T-matrix h as s everal desirable features. 
Up to corrections of order V no 3 , it reduces to the s-wave 
scattering length at T = 0. Thus the choice (fHH|l ~ lpHJ!)|) 
will reproduce the results for the weakly interacting Bose 



C. Renormalization Group Approach 



One very powerful method of quantum field theory 
is the Wilson renormalization group (RG) (Polchinski, 
ll984UWilson and Kogutl 11974ft . The basic idea is to sep- 
arate the momentum modes in the path integral into fast 
modes and slow modes by a cutoff. One then integrates 
out the fast modes. This yields an effective action for the 
slow modes, in which the coefficients of the original op- 
erators are renormalized and new operators are induced. 
By lowering the cutoff infinitesimally, one obtains a set of 
differential equations for the parameters in the effective 
action. Integrating out all modes down to k = yields 
the full effective action. 

Renormalization group techniques have been ap- 
plied to the homogeneo us Bose gas at finite tempera- 
ture by several au th ors {Andersen and Strickland! Il999t 
iBiilsma and Stooi Il996at iMetikas and Albert l2002|) . 
The first quantitative study was carried out by Bijlsma 
and Stoof. In that paper, they considered the one-loop 
diagrams that contribute to the effective chemical po- 
tential, effective four-point vertex etc. By introducing 
a cutoff, as explained above, a set of coupled differen- 
tial equations was derived. By solving these equations 
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numerically, they obtained the condensate density as a 
function of the temperature and the critical temperature 
as a function of the s-wave scattering length. They also 
derived equations for the fixed points and calculated crit- 
ical exponents. 

A somewhat different renormalization group 
approa ch was used by A n dersen and Strick- 
land. l)Andersen and Strickland. |l999). It is based 
on a derivative expansion of the effective action 
T[v] l)Morrisl 11994) . which is obtained by integrating out 
the quantum and thermal fluctuations. If one imposes 
an infrared cutoff k, one can expand the corresponding 
effective action [v] as 



r k [v] 



dr I d d x { V k (v) + L iZ £\ v )e i:j v~v 



a 



dr 



(196) 



where i,j = 1, 2, and repeated indices are summed over. 
Vk is the effective potential and zjp and Z^ are wave- 
function normalization constants, is the Levi-Civita 
symbol and Vi is the zth component of condensate v. The 
dots indicate all higher order terms in the derivative ex- 
pansion. Here and in the following, the subscript k indi- 
cates a dependence on the infrared cutoff. By lowering 
the cutoff k, one obtains a set of coupled integral equa- 
tions for the functions Vk, zj^\ Z^ ■ ■ ■. The leading or- 
der in the derivative expansion is defined by setting the 
coefficients Z^ and zjp to unity and the coefficients 
of all higher derivative operators in Eq. (|196|l to zero. 
This is called the local potential approximation. In the 
following, we restrict ourselves to the local potential ap- 
proximation and derive a flow equation for the effective 
potential Vfe(u). 



path integral into slow and fast modes separated by an 
infrared cutoff k. This is done by introducing a cutoff 
function Rk{p) or regulator, which we keep general for 
the moment. We add the term 

SkbPi,H = - f dT J dtx $R k (V-V*) 

x [V'iV 2 ^ + V2VV2] , (198) 

to the action Eq. I|121|) ■ The argument p of the function 
R k {p) has been replaced by y/— V 2 in coordinate space. 
After performing the shift (|37fl , the modified propagator 
reads 



D k (uJn,p) 



1 



p 2 (R k (p) + 1) + Vj 



p 2 (i? t (p) + i) + K ' + 2y V 



(199) 



where a prime on Vq denotes differentiation with respect 
to v 2 . The modified dispersion relation is 



e*(p) - ^/[p 2 (Rk( P ) + l) + V^ 

xiJ[p2(Rk(p) + l) + Vi + 2V>>v 2 



(200) 



By a judicious choice of R k (p), we can suppress the low 
momentum modes in the path integral and leave the high- 
momentum modes essentially unchanged. It is useful to 
introduce a blocking function fkip) which is defined by 



R k (p) 



1 - fk( P ) 
Ship) 



(201) 



1. One-loop effective potential 

We are now ready to calculate quantum and thermal 
corrections to the classical potential. We compute the 
one-loop effective potential which we will "RG improve" 
in the Section IIV.G.2I This method of deriving RG flow 
equations is conceptually and technically simpler than 
the direct application of exact or momen tum-shell RG 
techniques IjAndersen and Stricklandl 11999). 

The one-loop effective potential reads 

V = Vo + ^TtlogD- l (u n ,p) 

= - fiv 2 + lg V i + l^log[co 2 n + e 2 (p)} , (197) 

where the mean-field potential has been denoted by Vo = 
—[w 2 + 5<?i> 4 . We proceed by dividing the modes in the 



The blocking function satisfies 

lim/ fc (p)=0, lim/ fc (p) = l. (202) 

These properties ensure that the low-momentum modes 
are suppressed by making them very heavy and the high- 
momentum modes are left essentially unchanged. Typical 
blocking functions are shown in Fig. 1121 

The sharp cutoff function Rk{p) is defined by the 
blocking function fkip) = 9{p — k), and is shown in 
Fig. E| (solid line). It provides a sharp separation be- 
tween fast and slow modes. Using the sharp cutoff, the 
slow modes become completely suppressed in the path 
integral, while the fast modes are completely unaltered. 
The advantage of using a sharp cutoff is that certain inte- 
grals can be done analytically and the integro-differential 
RG-equation is reduced to a differential RG-equation. 

The solutions to approximate renormalization group 
equations (e.g. the truncation of the derivative expansion 
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the fast modes to the slow modes is completely ignored. 
Since all modes are integrated out independently, this 
is so metimes called the indepen dent mode approxima- 
tion (L iao and Strickland ll99E|) . The lack of feedback 
leads to a poor tracking of the effective degrees of free- 
dom. The situation is remedied by applying the renor- 
malization group, which effectively sums up important 
classes of diagrams. 



2. Renormalization group improvement 



FIG. 12 Sharp blocking function (solid line) and a typical 
smooth blocking function (dashed line). 



at some finite order) depends on the regulator function 
R k (p). There have been seve ral papers on the optimal 
choice of the cutoff fun ction ijAndersen and Strickland 
Il999t iLitiml 1200111200(1) . Here we consider the class of 
smooth blocking functions 



fJT(p) 



p" 



p 2 + k" 



(203) 



In the limit m — > oo, we recover the sharp blocking func- 
tion. In Scc. lIV.C.3l we shall see that a smooth regulator 
is better than the sharp one, but it comes at the expense 
of more complicated numerics. 

We return to the one-loop effective potential 
in l|197|l . Using the inverse propagator D7 (u) n ,p), the 
modified one-loop effective potential becomes 



V k = Vn 



(204) 



Upon summation over the Matsubara frequencies, we ob- 
tain 



V k = Vn 



d d p (1 , , 



(2tt 



I _ e -0ek(p) 



(205) 



The first term in the brackets is the T = part and 
represents the zero-point fluctuations. The second term 
includes thermal effects. Differentiation with respect to 
the infrared cutoff k yields: 



4* 



= -k 



d d p (dR k {p) 



(2tt) 



dk 



1 



x [p 2 (R k (p) + l) + V^ + V^ 



2e k {p) 



VI 



[l + 2n(e k (p))} 
(206) 



Eq. (|2U(j[l is an integro-differential equation for the one- 
loop effective potential. It is obtained by integrating 
out each mode independently, where the feedback from 



The easiest way of deriving the RG-improved ver- 
sion of Eq. (|2U6[1 is simply to make it self-consistent by 
replacing everywhere Vq by V k . This yields: 



k lk V * = " 



d d p (dR k {p)\ 1 



(2n) 



dk 



2e k (p) 

'VI 



[1 + 2n(e fc (p))] 



x [p 2 (R k (p) + l) + Vl + V^v'\ , (207! 
where the RG-improved dispersion relation is 



e k (p) = ^j[p 2 {R k (p) + l) + Vl] 

x^/[p 2 (R k (p) + l) + V k ' + 2V>>v 2 } , 



(208) 



and the primes in Eqs. I|207|l and l|208|l again denote 
differentiation with respect to v 2 . The self-consistent 
Eq. (|207|l is not a perturbative approximation, but is ex- 
act to leading order in the derivative expansion. This 
equation has been rigorously derived using the path in- 
tegral representa tion and the derivative expansio n of the 
effective action. ()Andersen and Stricklandl.ll999|) . There 
are other ways of regulating the one-loop expression( H97l) 
and subsequently "RG-improve" it. However, they do 
not always resum perturbation theory correctly (with 
all higher order diagrams included and all combinator- 
ical coefficients correct). These is sues have been studied 
in de tail by Litim and Pawlowski l)Litim and Pawlowskil 
I2002D . 

Note that since V k _ vanishes at the minimum of the 
effective potential, the dispersion relation in the broken 
phase reduces to 



£fco(p) = P^P 2 + 2V^ v 2 



(209) 



Thus, the Goldstone theorem is automatically satisfied 
for temperatures below T c . 

In the calculations that follow, we will restrict our- 
selves to using the sharp cutoff function. The integral 
over p in Eq. H207|l can be done analytically, resulting 
in a differential RG-equation. In this case, Eq. I|207|) 
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reduces to 



k-^V k = -\Sdk d \ e(fc) + 2Tlog 



1 - e 



-0e(k) 



> (210) 



where 



S d 



(2ir) d ' 



(211) 



and fid is the area of a (i-dimcnsional sphere whose ex- 
pression is given in the appendix. 

Since the factor dRk(p) /dk forces p — k, we have 
defined e(k) = Chip — k), where the dispersion relation 
e(k) is 
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(212) 



In order to solve H210|) , one must impose the correct 
boundary condition on the effective potential For k = 
co, no modes have been integrated out and Vk reduces to 
the classical potential Vq. In practise, one must impose 
the boundary condition at a large but finite value k — A. 

Eq. (|210|l has been solved numerically for the effec- 
tive potential Vk=o{v) with the above boundary condition 
in d = 3 dimensions for different values of T are shown in 
Fig. EH We have normalized the condensate v as well as 
the effective potential by the appropriate powers of the 
ultraviolet cutoff A. The curves clearly show that the 
phase transition is second order. For T < T c , the effec- 
tive potential has a small imaginary part, and we have 
shown only the real part in Fig. 1 1 31 The imaginary part of 
the effective potential does, however, vanish for T > T c . 
The effective chemical potential fik=o as well as the effec- 
tive quartic coupling constant gk=o are shown in Fig. 1771 
and both quantities vanish at the critical point. The 
corresponding operators are relevant and must therefore 
vanish at T c , and we see that the renormalization group 
approach correctly describes the behavior near criticality. 
Moreover, it can be shown that the sextic coupling g k 6 l 
goes to a nonzero constant at the critical temperature. 



3. Critical behavior and critical exponents 

In order to investigate the critical behavior near 
fixed points, we write the flow equation in dimcnsionlcss 
form using the dimcnsionless quantities 



(3 = 13k 2 , 
Vk = i3k~ d V k 



v = f3 1/2 k— v 
e(k) = k~ 2 e(k) 



(213) 



FIG. 13 Real part of the RG-improved effective potential 
I4=o (v) for various values of the temperature. The phase 
transition is clearly second order. 
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FIG. 14 The effective chemical potential /ik=o and the ef- 
fective quartic coupling gt=o near the critical temperature. 
Both vanish at T c . 



Eq. I|210(l can then be written as 



% + ^S d Pe(k) 



-S d log 1 — e 



-/3e(fc) 



(214) 



The critical potential is obtained by setting to zero the 
derivative with respect to k in Eq. (|214|> . Expanding in 
powers of 0e{k), we obtain 



■- 2 (d-2)v^_ + d 



V k - --S d $e(k) 

-S d \og [0e(k)] 



(215) 



Taking the limit (3 — > and ignoring the term which is 



29 



independent of v leads to 



2 ov 



V k = - 



\s d 



log [1 



log [l + V^ + 2Vl'v 2 ] 



(216) 



This is the same equation as obtained by Morris ijMorrisl 
1998) for a relativistic 0(2)-symmetric scalar theory in d 
dimensions to leading order in the derivative expansion. 
Therefore, the results for the critical behavior at leading 
order in the derivative expansion will be the same as 
those obtained in the d-dimensional 0(2)-model at zero 
temperature. 

The above also demonstrates that the system be- 
haves as a d-dimensional one as the temperature becomes 
much higher than any other scale in the problem Thus 
the system goes from being d+ 1-dimcnsional at low tem- 
perature to being effectively d-dimensional at high tem- 
perature and this is often called dimensional crossover. It 
is also referred to as dimensional reduction. The nonzero 
Matsubara modes decouple and the system can be de- 
scribed in terms of a cl assical field theory for the n — 
modes in d dimensions ijLandsmarl [19891 . We return to 
this subject in Sec. where we discuss the calculation 
of the critical temperature of a dilute Bose gas. 

The RG-equation (|207(l satisfied by V k [v] is highly 
nonlinear and a direct measurement of the critical ex- 
ponents from the numerical solutions is very time- 
consuming. This becomes even worse as one goes to 
higher orders in the derivative expansion and so it is 
important to have an additional reliable approximation 
scheme for calculating critical expone nts. In the fol l owinj 
we perform a polynomial expansion ijNicoll et all Il97- 
of the effective potential, expand around v = 0, and trun- 
cate at Nth order: 



2 9kV 



N (2n) 

9 k 

n— 3 



n 



-v 
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The polynomial expansion turns the partial differential 
equation 12101) into a set of coupled ordinary differential 
equations. In order to demonstrate the procedure we 
will show how the fixed points and critical exponents are 
calculated at the lowest nontrivial order of truncation 
(N = 2). We write the equations in dimensionless form 
using Eqs. I|213j) and 
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k Vfc 
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We then obtain the following set of equations: 



k 8k fLk = - 2 f I k + S d Pg k [l + 2n(e(k))} (219) 



k ~dk 9k = ^ ~ ^ 9k + Sd @ 9k 



1 



2(1 -fi k ) 
/3n(e(fc))[l + n(e(fc))] 



[1 + 2n(e(fc))] 
(220) 



A similar s et of equations was first obt ained in by Bijlsma 
and Stoof ( Biilsma and Stool Il99 6a) by considering the 
one-loop diagrams that contribute to the running of the 
effective chemical potential, the effective quartic coupling 
constant etc. 



The equations for the fixed points are 



k-Jik = 



k dk 9k = ° 



(221) 



If we introduce the variables r and s through the relations 



/'-A: 



9k 



1 



A*fc 



1 



MfcJ 



(222) 



and expand the equations in powers of /3(1 — the 
RG-equations can be written as 



Or 

dk 
ds 
dk 



= -2 [1 + r] [r - S d s] , 
= -s[4-d + ir-9S d s] 



(223) 
(224) 



We have the trivial Gaussian fixed point (r, s) — 
(0, 0) as well as the infinite temperature Gaus- 
sian fixed point (—1,0). Finally, for d < 
4 there is the infrared Wilson-Fisher fixed p oint 
((4 - d)/5, (4 - d)/ (55 d )) llWilson and Kogutl 119741) . 



3 and linearizing Eq. (|223|l around 
we find the eigenvalues (Ai,Aa) = 



Setting d — 
the fixed point, 
(—1.278, 1.878). The critical exponent v is given by the 
inverse of the largest eigenvalue: v = 1/A 2 = 0.532. This 
procedure can now be repeated including a larger num- 
ber, TV, of terms in the expansion Eq. I|217|l . The re- 
sult for v is plotted in Fig. ^| as a function of the num- 
ber of terms in the polynomial expansion (dashed line). 
Our result agrees with that of Morris, who considered 
the r elativist i c Q(2 )-model in d = 3 at zero tempera- 
ture ijMorrisl Il998'). The critical exponent v oscillates 
around the average value 0.73. The value of v never 
converges as N — > oo, but continues to fluctuate. For 
comparison, we have also showed the results using the 
smooth regulator function (|203|l with m — 5. Clearly, 
the convergence properties has improved signficantly. 

Our result should be compared to experiment on 
4 He and the e-expans ion which both give a value of 
0.67 l)Zinn-Justinl ll989). The result v = 0.73 is a leading- 
order result in the derivative expansion of the effective 
action i|196|) . The next-to-leading order in the deriva- 
tive expansion involves a set of coupled equations for the 
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potential V& and the wavefunction normalization terms 
zj 1 ^ and Z^j\ T hese equations were derived by Mor- 



ris l)Morrisl 119981) . A calculation of the critical exponent 
gives a value of v = 0.65. One expects the critical ex- 
ponent v to converge towards 0.67 as one includes more 
terms in the derivative expansion. 
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FIG. 15 The critical exponent v as a function of number of 
terms N in the polynomial expansion. 

We close this Sec. by comparing the two ap- 
proaches discussed here ll Andersen and Strickland 1 199(1 
iBiilsma and Stool , Il996a(l . The renormalization group 
equation (|210(l that we obtained in the local potential 
approximation is highly nonlinear. We made the fur- 
ther approximation by expanding the effective poten- 
tial in a power series, which turns the partial differen- 
tial equation (|210|l into a set of coupled ordinary dif- 
ferential equations. These flo w equations were first ob- 
tained by Bijlsma and Stoof l|Biilsma and Stool [l996a) 
by considering the one-loop diagrams that contribute to 
the various n-point functions, while neglecting their mo- 
mentum dependence. However, the two sets of RG- 
equations differ if one includes the momentum depen- 
dence of the one-loop gr aphs and go beyond the lo c al po- 
tential approximation. fAndersen and Strickland Il999t 
iLiao and Strickland 1 1995[) . 



D. Beliaev-Popov Approximation 

In Sec. Illl.BI we discussed in some detail the Bcli- 
aev approximation at zero temperature. We recall that 
this approximation is defined by all one-loop diagrams 
contributing to the self-energies. It was shown t hat the 
resul ting dispersion relation is gapless. Popov l)Popovt 
1983) generalized the Beliaev approximation to finite 
temperature, and hence it i s often called the Bel iaev- 
Popov (B-P) approximation l)Shi and Griffinl . |l998|) . He 
calculated the one-loop diagrams in the limit of zero 
external energy and momentum. However, one has to 
be careful because the limits u) — ► 0, p — > and 
p — ► 0, u> — ► do not commute at finite temperature. 
Later, Shi and Griffin have given a detailed discussion of 
the B-P approximation and explicit formal expressions 
for the normal and anomalous self-energies at finite tem- 



perature and arbitrary external energy u and momen- 
tum p. An analysis of the B-P approximation similar to 
the one presented at zero temperature in Sec. IHII would 
probably be significantly simpler. For instance it fol- 
lows immediately from the finite temperature version of 
Eq. i|105|> that the chemical potential is free of the in- 
frared divergences that plague separately the expressions 
for En(0,0) and £12(0,0). The expression for the chem- 
ical potential is 



2 1 

1" = rjv + 7.9 



Ap 2 — 4/i + 6gv 
u; 2 +e 2 (p) 



gv 



-g [3/i,i + J-i ,-1] 



'2^9 I dp 



p 2 (2p 2 -2n + 3gv 2 



-n{e(p)) . (225) 



In the high-temperature limit, where Pgv 2 <C 1, one finds 



gv 



3 T 
87r 2 gv 2 



V / 2g 3 v 2 



(226) 



Thus at high temperature, the dimensionless expansion 
parameter has an extra factor of T/gv 2 compared to the 
zero-temperature case, where it is y g 3 v 2 . 

So far our description of the dilute Bose gas has been 
in terms of spontaneous breaking of the U(l) symmetry 
and the use of a chemical potential to ensure that the 
mean number of particles is constant. There exist alter- 
native approaches that are number-conserving. Such an 
approach was first i ntroduced by Giradeau and Ar nowitt 
as early as in 1959 l|Giradeau and Arnowittl 11959) . This 
approach is variational in nature and, like the HBF ap- 
proximation, it exhibits a gap in the disp ersion relation . 
However, later it was shown by Takano ijTakanol Il96l|) 
that the gap is removed if one takes into account cubic 
terms in the Ham iltonian i n a co nsistent manner. More 
recently, Morgan ijMorganl Il999fl has discussed in detail 
number-conserving approaches and generalized them to 
trapped Bose gases. The starting point is the Hamilto- 
nian including interaction terms describing binary col- 
lisions written in terms of pair operators that conserve 
particle number. These operators are defined by 



(227) 



where aj are the standard annihilation operators and 
00 = (1 + ajao) -1 / 2 ^. The Hamiltonian can then be 
divided into terms containing zero, one, two, three, and 
four creation/annihilation operators. The terms that 
contain up to two operators are then diagonalizcd, while 
the cubic and quartic terms are treated as perturbations 
using first and second order perturbation theory. The 
validit y of perturbation theory is the usual requirement 
that V na 3 <§; 1 at T — 0. At high temperature, it follows 
from l|226|l that the requirement be T y/ g 3 v 2 / gv 2 <C 1. 
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In the zero temperature limit, the calculations reproduce 
several known results presented in Sec. IIIII The Hamil- 
tonian can be used to calculate the ground state energy 
density £ (rather than the free energy density J-) and a 
calculation involving the quadratic terms yields the re- 
sult of Bogoliubov given in Eq. (|92|l . A second-orde r cal- 
culat ion reproduces the na 3 log(na 3 )-term of Wu (jWiJ . 
1959), but one is faced with a logarithmic ultraviolet di- 
vergent term which is proportional to na 3 . As we argued 
in Sec. IIII.EI the correct way to treat this divergence 
is to absorb it in the coefficient of the operator (ip*ip) 3 , 
which represents 3 — > 3 scattering. The approach also re- 
produces Beliaev's results for the phonon spectrum 11071) 
and the damping rate 2|> . Fin ally, we mention that the 
modified Popov appro ximation ( Hut chinson et oil Il998t 
iProukakis et a?J . ll998ft has also been examined within the 
frame work of number-conserving approaches ijMorganl 
1999). 



V. CALCULATIONS OF T c 

The critical temperature for an ideal Bose gas is 
given by Eq. (|17|l . A natural question to ask is: what 
is the leading-order effect of a weak two-body interac- 
tion on the critical temperature of a homogeneous Bose 
gas? This question has been around for almost fifty years, 
but only very recently has the issue been s ettled. It has 
been discu ssed in detail by Baym et al. l|Bavm et all . 
1999.2001) and we follow to some extent their paper. 
Various approaches to the pro blem have als o been dis- 
cussed very recently by Haque ijHaauel [2003) . 

In the following, we assume that the interaction is re- 
pulsive, which corresponds to a positive scattering length 
a. One might think that effects of a repulsive interaction 
is to decrease the critical temperature of Bose gas. For 
instance, the superfluid transition in liquid 4 He, takes 
place at a lower temperature than that of an ideal gas 
of the same density. However, liquid 4 He is not weakly 
interacting and it turns out that the leading effects of 
interactions in the dilute Bose gas is to increase T c . 

The first paper in which a quantita tive prediction 
appe ars, is from 1957 by Lee and Yang ijLee and Yand . 
11957ft . In that paper, the authors predict that the crit- 
ical temperature increases compared to that of an ideal 
Bose gas, and that the increase is proportional to y/a. 
Later t he same authors pred icted that the shift is linear 
with a iftLee and Yand .L1958). This prediction was purely 
qualitative since neither sign nor magn itude were given . 
A cou ple of years later, Glassgold et al. (jdassgold et a,l\ . 
1960) also predicted an increase of T c which is propor- 
tional to ^/a. A coup le of decades la ter, the problem was 
revisited by Toyoda ijTovodal Il982ft . He predicted a de- 
crease of the critical temperature which is proportional to 
y/a. Since the sign agrees with the measurements on 4 Hc, 
there seemed to be at least qualitative agreement between 



theory and experiments. Long ago, H uang claimed an in- 
crease in T c that is proportional a 3 / 2 ijHuand . Il964ft , and 
very recently he clai ms in another p aper that T c increases 
proportional to ^[a ijHuand . Il999ft . From this selection 
of papers, it is clear that there has been considerable 
amount of confusion about how T c depends parametri- 
cally on the scattering length. Common to these r esults 
as well other attem pts to calculate T c ijSchakel 1 19941 
IWilkens et all l2000(l . is that they are based on perturba- 
tion theory. However, Bose condensation in a dilute gas 
is governed by long-distance physics that is inherently 
nonperturbative. These issues will be discussed later. 

There have also been other approaches to the calcu- 
lation of T c . Stoof and Bijlsma have carried out renor- 
malization group calculati ons of the critical temperature 
l|Biilsma and Stooill996a|) . This approach was discussed 
in the previous section and it predicts that the leading 
shift is proportional to a log a. Since this approach is 
purely numerical, it is difficult to take the limit a — > 
and thus obtain the correct dependence on a in the di- 
lute limit. Very recently, a calculation of T c based on 
the exact renormal ization group by Led owski, Hassel- 
mann and Kopietz l|Ledowski et all 12003ft . They calcu- 
lated the momentum-dependent two-point function and 
showed that the leading behavior is proportional to the 
scattering length a. 

The first Monte Carlo simulations for hard-sphere 
bos ons in the low-densi ty regime were done by Griiter et 
al. ijGriiter et a/.l .ll997). In that paper, the authors pre- 
dict a positive linear shift after extrapolating to the limit 
a — > 0. This result was somewhat surprising since some 
early Monte Carlo simulations as well as the experiments 
on 4 He, show a decrease in T c due to repulsive interac- 
tions. Later, it has been shown rigorously using effective 
field theory methods tha t the parametric depend ence of 
T c indeed is linear in a i|Bavm et all ll999.200i)K Thus 
in the dilute limit, we can write 



AT C 



= cn 1/3 a 



(228) 



where c is a constant that is to be determined. The prob- 
lem of determining the constant c has been attacked by 
analytical as well as numerical methods in recent years. 
These methods include high-precision Monte Carlo sim- 
ulations, 1/iV expansions, self-consistens calculations in- 
volving summation of bubble and ladder diagrams, and 
variational perturbation theory. We will discuss these in 
Secs. lV.B.lHV.B.3l 



A. Hartree-Fock Approximation and Breakdown of 
Perturbation Theory 



In this subsection, we will briefly discuss the 
Hartree-Fock approximation and show that it predicts 
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no shift in the critical temperature. We approach the 
phase transition from above, so the condensate density v 
is zero. 

In an ideal gas, the number density of excited parti- 
cles is given by Eq. (|14fl . which can be written as 



n cx = ^3-33/2(2), 

Arp 



(229) 



where the function g(z) is the polylogarithmic function 



9l{z) 



00 r 



(230) 



and z — is the fugacity. 

The simplest way to include the effects of interac- 
tions is to include the self-energy in the Hartree-Fock ap- 
proximation. The Feynman graph is the tadpole diagram 
shown in Fig. 




interacting gases as ideal gases with modified parameters 
and thus predict the same T c . 

Calculations using the Hartree-Fock approximation 
have been carried out by Huang. He applies a virial ex- 
pansion to Eq. H232|) and obtains a change in the crit- 
ical temperature which is proportional to a 1 / 2 . How- 
ever, this is an artifact of the approximation, as can be 
seen by including more terms in the expansion. Thi s 
was discussed in some detail in IIBavm et alJ.ll 999.2001). 
A cor rect treatment is given in e.g. l|Pethick and Smithl . 
2002). Similarly, the sum mation of the ring dia- 
grams ijHaugset et all Il998l) in the effective potential 
does not change T c . The reason is that these diagrams 
are evaluated at zero external momentum and therefore 
merely corresponds to a redefinition of the chemical po- 
tential. 

We have seen that a leading perturbative calculation 
in the scattering length a gives no corrections to the crit- 
ical temperature of a dilute Bose gas. One might try to 
improve on this result by going to higher orders in per- 
turbation theory. The Feynman diagrams contributing 
to the self-energy at second order in perturbation theory 
are shown in Fig. 1171 



FIG. 16 Self-energy diagram in the Hartree-Fock approxima- 
tion. 

The tadpole is independent of the external momen- 
tum and the expression for the Hartree-Fock self-energy 
is 



Shf = 2g 



P - V 



+ (p 2 ~ M) 2 



2gn 



(231) 



where we have used Eqs. (|14H and (|A.30|I . Thus a particle 
with momentum p effectively has the energy e(p) = p 2 + 
2gn, where the term 2gn arises from the mean field of 
the other particles. We can now generalize Eq. I|229l) by 
writing 



^3/2 (e^" 2 - >) 



(232) 




FIG. 17 Two-loop self-energy diagrams. 

If we focus on the contribution from the n = Mat- 
subara mode, the diagrams read 



S*(0,p) 



-4g 2 T 2 



d d k 



d d q 

(27r) d J {2ir) d 



( fc 2- M )(g2_ M) 2 > 

d d k f d d q 



oh, 9 9 f d d k f d d q 



(k 2 - fi)(q 2 - [i) {\p + q + k\ 2 ~ (i) 



(233) 



(234) 



Eq. H232JI shows that we must increase the chemical po- 
tential by an amount A/i = Shf to keep the same number 
density as that of an ideal gas at the same temperature. 
It shows in particular that n must approach Ehf from 
below to obtain the critical number density at a given 
temperature. Thus the critical temperature remains the 
same. The conclusion is that including a constant mean- 
field shift in the single-particle energies cannot change 
the critical temperature of a Bose gas. This is an exam- 
ple of the fact that mean-field theories effectively treat 



where the chemical potential acts as an infrared cutoff 
in the integral. The left diagram is independent of the 
external momentum. If we use the mean-field criterion 
that fi — > at the transition, the integral is linearly diver- 
gent in the infrared. The right diagram depends on the 
external momentum k. For fi = 0, it is logarithmically 
divergent as the external momentum goes to zero. As 
one goes to higher orders in the perturbation expansion, 
the diagrams become increasingly infrared divergent for 
/i = 0. If we denote a generic self-energy diagram with n 
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loops by E n , we have i|Bavm et all Il999.200l|) 



. a 

At 



(235) 



This shows that perturbation theory breaks down in the 
critical region due to the infrared divergences. Physically, 
these infrared divergences are screened and this can be 
taken into account by summing certain classes of dia- 
grams from all orders of perturbation theory. Examples 
of this is summation of bubble or ladder diagrams. We 
return to this issue at the end of this section. 



B. Dimensional Reduction 



In Sec. rrvTH we saw that the renormalization-group 
equations at high temperature reduce to those of a three- 
dimensional 0(2)-symmetric theory. This is an example 
of dimensional reduction and for the dilute Bose gas, it 
can be understood as follows. In the imaginary time for- 
malism, the fields are decomposed into modes which are 
characterized by their Matsubara frequency u> n — lixnT . 
At distances much larger than the thermal wavelength 
and for temperatures sufficiently close to the critical tem- 
perature, the time derivative term for n ^ is much 
larger than both the kinetic energy term and the effec- 
tive chemical potential term. This implies that the non- 
static Matsubara modes decouple and the long-distance 
physics can be described in terms of an effective three- 
dimensional field theory for the n — mode. The fact 
that the long-distance physics associated with the phase 
transition is well separated from the typical momentum 
scale T associated with the nonzero Matsubara modes, 
makes effective field theory methods ideal to study the 
phase transition. 

The effective three-dimensional theory that de- 
scribes the long-distance physics can be constructed us - 
ing the methods of effective field theory l|Georgil Il993). 
Once the symmetries of the theory have been identified, 
one writes down the most general Lagrangian C e g that is 
consistent with these symmetries. In the present case, 
we simply have a complex scalar field with an 0(2)- 
symmetry. In addition, there a is three-dimensional rota- 
tional symmetry. The effective three-dimensional theory 
is then described by the action 



[V (</>*</>)] 2 and (<p*(f>) 3 . The relation between the param- 
eters in the effective theory and in the full theory can be 
determined by perturbative matching; one requires that 
the effective theory (|236|l reproduces static correlators at 
long distances R 3> 1/T to a specified accuracy. The rea- 
son why the parameters of the effective three-dimensional 
theory can be determined in perturbation theory, stems 
from the fact that the coefficients of the effective theory 
encode the short-distance physics at the scale T which is 
perturbative and that the matching procedure does not 
involve the nonpcrturbative long-distance physics. 

At the tree level, the matching can be done simply by 
inspection. By comparing the action H121[> that describes 
the full four-dimensional theory with the action l|23t)|) 
that describes the effective three-dimensional theory, we 
can read off the relation between the fields in the two 
theories. This yields 



(237) 



Similarly, by matching the other terms in the action in 
the two theories at the tree level and using l|237|) . one 
easily finds 



M3 

u 



3gT 



(238) 
(239) 



Higher order operators as well as corrections to the co- 
efficients of the operators in Eq. i|23fci|) can be ignored 
at the order of interest in the diluteness expansion. For 
instance, the coefficient of the operator [V (<fi* <fi)] 2 is pro- 
portional to a 2 At- From dimensional analysis, it follows 
that the contribution to physical quantities from this op- 
erator is then suppressed by a factor of n 1 / 3 a compared to 
the operator (^>*</>) 2 . The contribution to physical quan- 
tities from other operators are analyzed in a similar man- 
ner. 

For an ideal Bose gas, the critical temperature is 
given by Eq. 1171) . Equivalently, the critical number den- 
sity w? c at fixed temperature satisfies n°A^ = £ (§). Due 
to the repulsive interactions in the dilute Bose gas, the 
critical number densitiy changes. The first-order change 
in the critical temperature AT C = T c — T® is related to the 
first-order change An r = n a — njj in the critica l number 
density at fixed T c by l|Bavm et all 11999 200l|k 



Soft = 



d 3 x 



24 



, (236) 



AT C 2 [n c (T c ) - n° c {T c )] 

T c ° ~~ 3 n 
1 T°A(</)*<j)) 
~3 n° ' 



(240) 



where we have used the conventional normalization of 
an 0(2)-invariant theory. The dots indicate operators 
with more derivatives and more fields. Examples are 



The factor 2/3 in the first line comes from the relation 



(n 1 



Q\2/3 



The last equality follows from the fact 



that n = (il , *tp) an d that the contribution from the the 
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zeroth Matsubara mode is hT((j)*(p), which follows from 
Eq. 123711 . 

In order to calculate the critical temperature, we 
must evaluate the quantity A((f)*(p). We discuss this next. 



1. 1/N expansion 



FIG. 18 Three-loop Feynman diagrams contributing to 
A(0 2 }. 



The 1/N expansion is a nonpcrturbative method 
thas has been widely us ed in high-energ y and 
condensed matter physics zin and Wadia, 1993; 
Moshc and Zinn- Justin! l2003|) . In condensed matter 
physics, it has been used to study the critical behavior 
of O( N) spin mod e ls and calculate their critical expo- 
nents l)Zinn- Justinl GUI. The idea is to generalize a 
Lagrangian with a fixed number of fields to N fields and 
then let N be a variable. The expansion is defined as 
an expansion in powers of 1/N while gN is held fixed (g 
is the coupling constant). The method is nonperturba- 
tive in the sense that calculations at every order in 1/N, 
there are Feynman diagrams contributing from all orders 
of perturbation theory. In this way, one sums up graphs 
from all orders of perturbation theory. One hopes that 
this expansion captures some of the essential physics that 
cannot be captured by e.g. perturbative methods. 

The critical temperature was recently calculated by 
Baym, Blaizot, and Zinn-Justin ijBavm et [2000) m 
the large- N limit. The next-to-le ading-order result was 
obtai ned by Arnold and Tomasik l|Arnold and TomasikL 
120001) . 

In the present case, the Lagrangian (|236|l is gener- 
alized to a scalar field theory with N real components. 
The Lagranigan is now O(N) invariant and reads 



1 ^2 1 , 

2^ + 24 U( 



(241) 



where i runs from 1 to N. Summation over i is implic- 
itly understood. The large- N limit is obtained by taking 
N — > oo, while keeping uN constant. 

Generally, the diagrams that contribute to A(0 2 ) 
can be obtained from vacuum diagrams by inserting an 
operator <f) 2 . For example, at the three-loop level, there 
are two diagram that contribute and they are shown in 
Fig diagrams in Fig. ^] It can be shown that the first 
diagram is suppressed by a factor of 1 /N relative to the 
second. 

The Feynman diagrams that contribute to A(</> 2 ) at 
leading order in 1/N are shown in Fig. ^5] The dot de- 
notes an insertion of the operator <f> 2 . These diagrams 
are called bubble diagrams and the bubble summation is 
thus exact in the large- N limit. The expression for the 




FIG. 19 Feynman diagrams contributing to A((f> 2 } to leading 
order in 1/N. 



expectation value A(</> 2 ) then becomes 



-N 



d d p 1 
(2tt)<V 



S(p) , 



(242) 



where £(p) is the self-energy. The Feynman diagrams 
for the self-energy are obtained from those in Fig. I19l bv 
cutting the propagator line that goes through the dots. 
The self-energy diagrams are shown in Fig. 1201 




FIG. 20 Feynman graphs contributing to the self-energy to 
leading order in 1/JV. 



After mass renormalizatio m so that S(0) = , the 
expression for the self-energy is ijZinn- Justinl . Il989|) 



e(p) 



d d k 



N J (27r) d 6/Nu + T^k) 



| P + k| 2 k 2 



(243) 



where T\ (k) is the one-loop contribution to the four-point 
function: 



Tt{k) 



d d q 



1 



(27r) d <? 2 |q + k| ; 



(244) 



In the appendix, we show how to calculate the function 
T\(k) in dimensional rcgularization. One finds 



Tx{k) = M 2 



l)r(| 



ud-4 



(245) 



2 2d- 37r ^ r (<^!) 
The next step is to evaluate the integral over p in 
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Eq. (212)1 . Using Eq. (A~2"T)) . we obtain 



d d p 1 
(27r) d p 4 



|p + k| 2 /c 2 



M 



2c 



r( 3 -|)r(|-2) , d _ 



(246) 



Inserting Eqs. I|245(l and Ij246|) into Eq. H243|) . we obtain 
T(3-f)r(f-2) 



-M 



d d k k d - 6 
(2ir) d a + bk d ~ A ' 



where 

a = 6/Nu , 
b = M 
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The final step consists of integrating over k. Using 
Eq. i|A.22|) in the Appendix we obtain 



d-2 3-d 



A(</> 2 ) = M 4e a^b^ 

r(3-f)r(f-2)r(|^)r(| E f) 



2 3rf - 6 vr d -^r(^) 



The limit d — > 3 is regular and one obtains 
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Inserting the result (|251|l into (|240|) , we obtain the critical 
temperature to leading order in l/N: 



AT C 
TP 



8tt 



3[C(1)] 
2.33n 1/3 a 



4/3 



1/3 



(252) 



where we have used that u = 3gT and g = 8ira and set 
N = 2. Thus T r i ncreas es linearly with a. As noted be- 
fore l)Bavm et all 12000). the result (|252|l is independent 
of N. However, it is only valid in the limit N — > oo. 

The l/N correction to the above result 
has rece ntly been calculat e d by Arnold and 
Tomasik (Arnold and Tomasikl 120001) . This is a 
very lengthy and technically complicated calculation. 
For instance the integrals that one encounters at order 
l/N are difficult to evaluate in 3 — 2e dimensions. 
Instead, Arnold and Tomasik always reduced their 
diagrams to unambigious integrals in three dimensions 



that are simpler to evaluate. We shall not review 
the calculation, but merely state the result. Through 
next-to-leading order in l/N, the shift in T c is 



AT C 
TP 



Hem 



3\l4/3 



1 - 



0.5272 



N 



(253) 



For N = 2, it gives a correction of only 26%: 
AT r _ „ 1/3 



TP 



Uln^a . 



(254) 



2. Monte Carlo simulations 



The action (|236|l can also be used as the starting 
point for numerical calculations of the critical temper- 
ature and other nonuniversal effects of a dilute Bose 
gas. This is done by putting the theory on a lattice and 
using Monte Carlo methods to solve the theory non- 
perturbatively. Such numerical si mulations have been 
carrie d out by sev era l gro up s ll Arnold and Moore , 
20011 I Arnold oil l2002t ICxriiter et all 11997 : 
Holzmaim^m^lCrauthL Il999t iKaghirrnikov et cd. , 
20011) . This approach has been discussed in considerable 



detail by Arnold and Moore (Arnold et a/1 l2002t) . The 
valu e of the coef fi cient in (|228|l reported by Gruter et 
al. (Gruter et all Il997f) is c ~ 0.34 while Holzmann 
and Krauth (Holzmann and KrauthL [l999) obtained 
c ~ 2.3. The most recent values reported by Arnold and 
Mo ore { Arnold and Mppre|| 200ll) and by Kashurnikov et 
al. (Kashurnlkov^^d 12001)) are c ~ 1.32 and c ~ 1.29, 
respectively. One source of discrepancy lies in the non- 
linear corrections to T c as a f unction of a at the de nsities 
where the simulations of llGriit er et all 1 19971) wer e 
carried performed. In (HolzmamHm'cH<xautbt Il999|) . 
the authors expand the integrand in the path integral 
in powers of the interaction and keep only the first 
term in that expansion. This perturbative treatment of 
the interaction is incorrect since the physics close the 
phase transition is inherently nonperturbative. This is 
in contrast to the Monte C arlo simulations based on 
the action H236H c a rried out (Arnold and Moorel l200ll 
iKashurnikov et al\ . 12001)) . These calculations agree 
within error bars. We will regard these lattice results as 
the correct result for the coefficient c in Eq. (|228|l . 



3. Other calculations 

In this subsection, we will briefly discuss other calcu- 
lations of the shift in T c in the dilute Bose gas. Generally, 
the expectation value (4> 2 } can be written as 



d d k 



1 



(2w) d fc 2 + (i 3 + E(fc) 



(255) 
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where E(fe) is the self-energy function. The critical point 
is determined by the condition that the correlation length 
becomes infinite, or equivalently that the effective chem- 
ical potential /13 + E(0) vanishes: 



M3 + S(0) = 



(256) 



In absence of interactions, the condition l|256(l reduces to 
the well known /13 = 0. Using Eq. (|256|) to eliminate the 
chemical potential (13, we can write Eq. I|24(J|) as 



d d k 



1 



1 



E(fc) - E(0) k 2 



(257) 



The next step is to make an approximation for the self- 
energy function appearing in Eq. (|257|l . In the previ- 
ous section, we used the la rge- TV expression for the self- 
energy Baym et al. (|Bavm et all 11999.200 1|) con- 
sidered three different equations for the self-energy. 

• One-bubble approximation: 

The self-energy is approximated by the second diagram 
in Fig. [T7| 



E(p) - E(0) = -2 9 2 T J *JL Tl 



(k) 



1 



|p + k| 2 + E(|p + k|) - E(0) k 2 + E(fe) - E(0)_ 

(258) 

1 Bubble-summation approximation: 



The self-energy is approximated by the bubble sum in 

Fig. na 



d d k Ti(fc) 



2gTt{k) 
1 1 



|p + k| 2 + E(|p + k|) - E(0) k 2 + E(fc) - E(0). 

(259) 

' Ladder-summation approximation: 



The self-energy is approximated by the ladder sum sim- 
ilar to the bubble sum. 

w ^ f ddk T i( fc ) 



1 



1 



|p + k| 2 + E(|p + k|) - E(0) k 2 + E(fc) - E(0)_ 

(260) 



Note that Eqs. (|258(l - l|260|) have been made self- 
consistent by replacing the free propagators on the right 
hand side by the interacting propagators. Note also that 
the only difference between the bubble sum and the lad- 
der sum is a factor of two. 

Eqs. (|258|I " H26U|) have been solved numerically and 
the results were used to evaluate Eq. H257|) to obtain the 
corresponding shifts in T c . The results for the coefficient 
c are 3.8, 2.5, and 1.6, respectively and thus within a 
factor three. Note also that the prediction for the shift 
in T c from the non self-consistent bubble summation (the 
leading 1/JV result with N = 2) is very close to the result 
from self-consistent bubble summation. 

We can gain more insight into the mechanism behind 
the increase in T c by looking at the modification of the 
spectrum for small momenta. Let us consider the non 
self-consistent bubble sum: 



;(fc) = fc 2 + E(fc)-E(0) , 



where 



m) = f ddk T *( fc ) 1 
1 ' J {2^1 + gT x {k) | P + 
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For small moment a k, the difference E(fc) — E(0) can be 
approximated by ijBavm et all Il999.200l)) 



E(fe) - E(0) 



2 , 2 /. k 



k < k c , (263) 



where k c = 8ir 2 a/X^ is the screening wave number. The 
logarithmic term in (|263|l indicates a modified spectrum 
for small wave numbers ~ k 2 ~ v and thus a harden- 
ing (The spectrum is of the form fc Q , where a < 2) of 
the spectrum compared to the noninteracting case. The 
hardening results from correlations among particles with 
low momentum, which leads to a decrease in the critical 
density and thus an increase in the critical temperature. 
Other approximations show a different functional depen- 
dence of the difference E(fc) — E(0) at small k, but the 
basic mechanism remains the same, namely a hardening 
of the spectr um at the critical temperatu re. It has been 
pointed out ijBavm and Holzmannll2003T) that the hard- 
ening of the spectrum for the n = Matsubara mode only 
takes place exactly at T c . The Bogoliubov operator in- 
equality guarantees that the spectrum remains quadratic 
away from T c . 



me thods have als o 
( Braaten and Radesctl 



been 



Variational 

used recently 

Ide Souza Cruz et all l200ll l2002j) to calculate T c . 
The basic idea is to compute (125711 using an effective 
three-dimensional Lagrangian that has been reorga- 
nized according to the discussion in Sec. IIV.B.2I The 
Lagrangian is written as C = Cq + C m t , where 



Co 



2 j,2 



-m 4> ; 



(264) 
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Ant = - 5(H3 



(265) 



and i runs from 1 to N. Calculations are carried out by 
using 5 as a formal expansion parameter, expanding to 
a given order in S and setting S = 1 at the end of the 
calculation. Finally, we need to give a prescription for 
the mass parameter m. In the calculations below, we 
will use the PMS criterion. In this context, it reads 



in three dimensions. The integral is calculated in the 
appendix. Using Eq. (|A.26|I . we obtain 



E 2 (fc) - Ea(0) 



(4tt) 2 



3m k 1 j k 2 + 9m 2 

k 3m 2 ^ 9m 2 



(270) 



0A< 



9m 



= 



(266) 



If we wish to apply variational methods, to calculating 
the shift in T c , we need to generalize the quantity A(</> 2 ) 
appearing in Eq. (|257|l to the field theory defined by 
Eqs. 1)2641) and (|265|) . One must be able to expand this 
quantity in powers of S and it must reduce to A(<j) 2 } when 
5=1. Several prescriptions for g eneralizing A(4> 2 ) have 
been proposed in the lite r aturel Br aaten and Radescil 
l2002t Ide Souza Cruz et oil 120011 120021) . Some of these 
prescriptions are well behaved in the limit N — > oo and 
some of them are not. Since the result for the shift in T c 
is known analytically in this limit, this a desireable prop- 
erty of a prescription. Two gener alizations that have this 
prope rty have been considered l)Braaten and Radescul 
2002): 



A h (cb 2 ) 



N 



N 



d d k 
{2ir) d 

d d k 
(2n) d 
1 



, (267) 



1 1 

k 2 + S(fc) - S(0) ~ k 2 

1 

k 2 + m 2 (l - 5) + S(fc) - S(0) 

(268) 



k 2 +m 2 (l-S) 



In the following, we will consider A a ((j) 2 ). The strategy 
is to calculate the difference £(fc) — 2(0) in a powers se- 
ries in S, substitute the result into Eq. (|267|l . and finally 
expand the resulting integral in powers of d. The Feyn- 
man diagram that contributes to the self-energy to first 
order in 8 is the leftmost diagram in Fig. El It is inde- 
pendent of the external momentum and so the difference 
S(fc) — S(0) vanishes. The first nonzero contribution to 
the quantity E(fc) — S(0) is then given by the two- loop 
diagram in Fig. 1171 The expression is 



e 9( o) = s 2N{N+2 ^- 2 r d3 p f ^ 



6 



1 



1 



(2tt) 3 J (2tt) 3 
1 



p 2 + m 2 q 2 + m 2 (p + q + k) 2 + m 2 
1 1 1 



p2 _|_ m 2 ^2 _|_ m 2 ^ _|_ _|_ 



, (269) 



The self-energy (|270H is itself second order in 6. The 



second-order result for A a 



is then obtained by ex- 



panding H267|) in powers of the subtracted self-energy and 
keeping only the first term: 



A^(6 2 ) = 



d 3 k 1 



[Ea(fc) - E 2 (0)] , 



(271) 



where the superscript indicates the order of S. Using 
Eq. JX271 



M 2 Hd> 2 ) = 



i 



i 



108(47r) 3 m 



S 2 u 2 N(N + 3) 



(272) 



The PMS criterion l|266|l has no solution at second order 
in 5 because (|272|l is a monotonic function of m. Thus 
one has to go the third order in order to obtain a value for 

m. T he result for A^ ((f) 2 ) reads (jBraaten and Radescul 
2002) 



Ai 3) <^> 
where 

h = - 

Here 



24(4tt) 



^— ^— S 2 u 2 N(N + 3)(l + -S 

108(4tt) 3 m V ; V 2 

^ 3 ^(^ + 2)(^ + 8)w 3 / 3 , 



7T 2 + 161og- + 12Li 2 (-l/3) 



(273) 



Li„(x) = 



i=i 



(274) 



(275) 



is the polylogarithmic function. At this order, the PMS 
criterion has a single real solution: 



1.04u 



A^ + 8 
24tt 



(276) 



The resulting value for A^ {4>1) IS 



where the subscript n indicates the order in the loop 
expansion. We have set d = 3 since the integral finite 



A^itf) = 0.4813 



iV+2 
iV + 8 



Nu \ 
96tt 2 J 



(277) 
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Setting N = 2, the prefactor bec omes 0.19. The lattice 
results l)Arnold and Moorel is 

A(tf) = (0.284 ±0.004) (~^) ■ (278) 

Thus the third-order approximation differs from the lat- 
tice M onte Carlo result by 66%. Similarly, it was 
found ([Braaten and Radescm I2002R that the fourth or- 
der approximation differs from the numerical simula- 
tion s by 61% In the paper b y Braaten and Rade- 
scu (|Braaten and RadescrJ . 12002') , the convergence of the 
linear 6 expansion to the exact result in the large- N limit 
was studied as well. It was shown that it converges to 
the lattice Monte-Carlo result l|252ll . but that the con- 
vergence is rather slow. 

A straightforward application of the linear delta- 
expansion to field theor etic problems have recently been 
criticized by Kleinert ijKleinertl l2003albft . O ne must 
take into account the correct Wegner exponent l)Wegneii 
Il972j) that governs the approach to the strong-coupling 
limit. The method is then called variational perturba- 
tion theory (VPT) l|KleinertL ll998T) . The failure to take 
into account the correct Wegner exponent is the reason 
for why one finds complex-valued solutions in the lin- 
ear delta-expansion. These matters are highly technical 
and beyond the scope of this paper. Interested read- 
ers are re ferred to the textbook by Kleinert an d Schulte- 
Frolinde l|Kleinert and Schulte-Frolindell200ll) . VPT has 
very reently b een applied in a seven-loop calculation 
by Kastening l lKaj^gnineL l2003bl) (See also ( Kas tenind . 
2003 a| iKleinertl l2003a|) for five and and six-loop calcula- 
tions.). The result for the coefficient is c = 1.28 ± 0.1 in 
excellent agreement with lattice field theory results. We 
also note that seven-loop calculations have been carried 
out for the case N = 1 and N — 4 as well. The values for 
the coefficient c are 1.07 ± 0.10 and 1.54 ± 0.11, respec- 
tively. These results are al so in good agreement with the 
lattice simulation of Sun ijSunl 2003) who obtained the 
values 1.09 and 1.59, respectively. 

We close this section by listing in Table. [I] the pre- 
dictions for the shift in the critical temperature that has 
been obtained by the various methods discussed in this 
section. 



VI. CONCLUSIONS 

In the present paper, we have extensively discussed 
the dilute Bose gas at zero and finite temperature. Using 
effective field theory methods, a systematic perturbative 
framework that can be used to calculate any quantity of 
the dilute Bose gas at zero temperature was set up. For 
instance, many of the classic results for the weakly inter- 
acting Bose gas are derived in an efficient and economi- 
cal manner. Nonuniversal effects are another application 
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= 2.33n 1/3 a, 


Leading order 1/N 
(rjavm et al.. zUUUi. 


AT C 
To 


— X. ( 171 (Ij 


Next-to- leading order 1/iV 
f Arnold and Tomasik. 2000). 


AT C 
To 


= (1.32 ± 0.02) n 1/3 a, 


Lattice (Arnold and Moore. 2001"). 


AT C 
To 
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Lattice f Kashurniko^^^^^/Z^. 2001 1 . 


AT C 
To 


= 0.7n 1/3 a, 


One-bubble approximation 
(Holzmann et al. 1999). 


AT C 
To 


= 3.8n 1/3 a, 


One-bubble self-consistent 
approximation 
I Ddvni et at., iyyy,zuuij. 


AT C 
To 


= 2.5n 1/3 a, 


Ladder- summation approximation 
( rSavm et at.. iyyy.zUUlj. 


AT C 
To 


1 fin 1 / 3 /! 


Bubble summation approximation 
fBavm et al. 1999.2001). 
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To 


= (1.23±)n 1/3 a, 


RG in three dimensions 
CLedowski et al. 2003). 


AT C 
To 


= (1.3±0.4)n 1/3 a, 


Simulations of classical field theory 
(Davis and Morean. 2003). 



TABLE I The critical temperature for a dilute Bose gas ob- 
tained by various analytic and numerical methods. 



where effective field theory methods are ideal. For in- 
stance, it was used to solve the long-standing problem of 
calculating the full order-na correction to the ground 
state energy density of a weakly interacting Bose gas. 
Similarly, it would be of interest to calculate higher-order 
corrections to the condensate density an d compare it with 
numerical results l|Giorgini et a/.lll999j) . The strength of 
the effective field theory approach lies in the fact that 
it is a systematic approach. To any given order in the 
low energy expansion, only a finite number of terms in 
the effective Lagrangian contribute to a physical quantity 
and this allows us to make definite predictions. Although 
effective field theories are nonrenormalizable in the tra- 
ditional sense of the word, one can carry out renormal- 
ization systematically order by order in the low-energy 
expansion. 
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We have also discussed different approaches to the 
dilute Bose gas at finite temperature. Although both the 
Bogoliubov and Popov approximations are gapless, they 
break down at sufficiently high temperature. Since the 
interactions between excited bosons are ignored in the 
Bogoliubov approximation, it is only valid at very low 
temperatures, where one can neglect the thermal deple- 
tion of the condensate. The Popov approximation can 
be used at much higher temperatures. However, being a 
mean-field theory, it breaks down in the critical region. 
It incorrectly predicts a first-order phase transition with 
the same T c as the ideal Bose gas. The full HFB approxi- 
mation violates the Hugcnholz-Pincs theorem, which any 
reasonable approximation should incorporate. The mod- 
ified Popov approximation based on the many-body T- 
matrix is an improvement over the original Popov ap- 
proximation because it correctly predicts a second-order 
phase transition. The r enormalization group equations 
that have been derived l|Andersen and Strickland! . Il999t 
iBiilsma and Stool Il996a|) show that the critical proper- 
ties of the dilute Bose gas are give by a three-dimensional 
spin model with a continuous 0(2) symmetry. Explicit 
numerical calculations demonstrate that the system un- 
dergoes a second order phase transition where the ef- 
fective coupling constant vanishes at the critical point. 
Of the other approaches to the thermodynamics of the 
dilute Bose, optimized perturbation theory is probably 
the most promising. It is a systematically improvable 
expansion with significant flexibility with respect to the 
choice of parameters. This approximation also respects 
the Goldstone theorem order by order in the perturbative 
expansion. To go beyond one loop for the free energy in 
optimized perturbation theory is a difficult, but not im- 
possible task. The complicated setting sun diagrams that 
appear at the two-loop level arc most easily calculated 
using the method where one separates the diagrams into 
contributions from z ero, one, and two Bose-Einstein dis- 
tribu tion functions ijAndersen. Braaten and Strickland! . 
1200(1) . 

The problem of calculating the shift in the critical 
temperature of a weakly interacting Bose gas has been a 
long-standing problem. Using effective field theory meth- 
ods to obtain an effective three-dimensional field theory 
combined with high precision lattice calculations has set- 
tled the issue in a very elegant way. The next-to-leading 
order result in the 1/7V expansion, shows that this ex- 
pansion works surprisingly well for T c ; for N = 2, the 
result is 29% higher than the predictions from the lattice 
simulations. Among other notable calculations, we men- 
tion the application of variational perturbation theory 
through seven loops. The agreement with lattice data 
for N = 1, 2 and 4 is convincing. 

This summarizes our current understanding of some 
aspects of the homogeneous Bose gas at zero and finite 
temperature. Significant progress has been made in the 
last ten years, but there are still open problems. Wc 
hope that this review has stimulated the reader to further 



research in the field. 
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APPENDIX A: Calculational Details, Notation and 
Conventions 



In this appendix, we give some calculational details 
that may be useful for the reader who is interested in 
going through the calculations in detail. We also define 
our notation and conventions used throughout the paper. 



1. Zero temperature 

All the zero-temperature calculations are carried out 
in real time. Loop integrals are integrals over real en- 
ergies u> and over three-dimensional momenta k. The 
integrals over u> are evaluated using contour integration. 



The specific integrals needed are 



^-lo g [^-e 2 {p)] 



duj 



1 



i 



2tt [tu 2 -e 2 {p)} 
du 1 



2e(p) 
i 



2ir [uj 2 -e 2 (p)} 2 4e 3 (p) 



(A.l) 
(A.2) 
(A.3) 



Some momentum integrals are divergent in the in- 
frared or in the ultraviolet, or both. Dimensional regu- 
larization can be used to regularize both the ultraviolet 
and infrared divergences in three-dimensional integrals 
over momenta. The spatial dimension is generalized to 
d = 3 — 2e dimensions. The continuum limit is taken by 
replacing sums over wave vectors by integrals in d — 3— 2e 
dimensions: 



V ^ 
p 



M 



2 c 



d d P 



(A.4) 



Here, M is a renormalization scale that ensures that the 
integral has the canonical dimension also for d ^= 3. In the 
following, we absorb the factor M 2e in the measure and so 
it will not appear explicitly. Integrals are evaluated at a 
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value of d, for which they converge and then analytically 
continued to d = 3. 



The integral I m , n is defined by 
d d p 



l (A)=M i 



P 



(2w) d p n (p 2 +A 2 ) n / 2 



(A.5) 



With dimensional regularization, I m ,n is given by the for- 
mula 



^771,71 (A) 



^ d j^j2e a d+2m-2n 



(2ir) d 



I) 



2T(f) 



(A.6) 



where f2 a 
sphere. 



2 7r d / 2 /r (|) is the area of a d-dimensional 



The integrals I m , n satisfy the relations 
dA 2 — 2 m ^i ,n + 2 ' 



(e? + 2m — n) I n 
A 2 I 



nl. 



m+2,n+2 j 



= I, 



m — 1 ,n— 2 m + 1 , n 



(A.7) 
(A.8) 
(A.9) 



The first relation follows directly from the definition of 
Im,n- The second relation follows from integration by 
parts, while the last is simply an algebraic relation. 

The specific integrals we need in the calculations are 
listed below. In the limit (i — s- 3, they become 



Io.-i 
h,i 

h, 5 
h,3 
Io,i 

1-2,-1 
1-1,1 



A 5 



157T 2 ' 
A^_ 

3tt 2 ' 

_ A^_ 

6tt 2 ' 
4A 

A 

A 

2^2 ' 
A 

1 

~4A? 



4 - L - 7 + log(7r) 



2 - L - 7 + log(7r) 



(A.10) 
(ATI) 
(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 

(A.17) 



where L = log(A 2 /M 2 ). The integrals Z-i,i and J_2, -i 
are both logarithmically divergent in the infrared and 
this shows up as a pole in e. 

We also need to calculate some integrals in d dimen- 
sions that depend on the external momentum k. The 



integral Ti(fc) defined in Eq. I|244[) is 
d d q 1 



(27r) d 9 2 |q + k| 2 



(A.18) 



By introducing a Feynman parameter y, we can write the 
integral as 



T(k) 



d d q 



dy J (2n) d ( g 2 +m 2 ) 2 



(A.19) 



where m = y{\ — y)k 2 . First integrating over q, and 
then over j/ gives 



TO) = Af 2 



T(2-|) f 1 
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m 
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(A.20) 



Other integrals needed are calculated in the same man- 
ner. Specifically, we need the integrals 



d d p 1 
(2Tr) d p^ 
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|p + k| 2 k 
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1 \ 4-cZ i 1 1 4-d ) -2 3-d 



-M 



2 d -%*r(f) 



(A.22) 



Note that the integral in (|A.21ll vanishes in d = 3 dimen- 
sions due to the factor T(— e) in the denominator. The 
integral (|A.22|) has a pole in e. When these integrals are 
combined, the limit d — > 3 is regular. 

We need to evaluate the subtracted two-loop self- 
energy in the linear delta expansion. The expression is 

d d p f d d q 



1 



1 



1 



p2 _|_ m 2 g2 _|_ m 2 _|_ g _|_ 7^2 _|_ m 2 

1 1 1 



p 2 + 77i 2 g 2 + m 2 (p + q) 2 + m 2 



(A.23) 



This integral is ultraviolet finite and can be calculated 
directly in three dimensions by going to coordinate space: 

Ea(ft) - S 2 (0) = J d 3 r [e lkr - l] V 3 {r) , (A.24) 
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where V(r) is the Fourier transform of the propagator: 



V(r) 



d d P 



(27r) d p 2 + m 2 



(A.25) 



This is the usual Yukawa potential. Integrating over r 
yields 



E 2 (fc)-£ 2 (0) 



1 



(47T) 



3m k 

1 — arctan 

k 3m 



1 k 2 + 9m 21 

2 ° g 9m 2 
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Finally, we need to evaluate the integral 
d d k 1 



A( 2 )(</> 2 ) = 



(27r) d /c 4 



[S 2 (fc) - E 2 (0)] 



The integral is finite in the ultraviolet and can thus be 
evaluated directly in three dimensions. One obtains 



Ai 2 V> = - 



1 1 

6(47r) 3 m 



(A.27) 



2. Finite temperature 

In the imaginary-time formalism for thermal field 
theory, the 4-momentum P = (w n ,p) is Euclidean with 
P 2 = w 2 + p 2 . The Euclidean energy po has discrete 
values: uj n — 2nirT for bosons, where n is an integer. 
Loop diagrams involve sums over u n and integrals over 
p. With dimensional regularization, the integral is gen- 
eralized to d = 3 — 2e spatial dimensions. We define the 
dimensionally regularized sum-integral by 



= M 2e T 



d d p 



l= 2ri7rT J v ; 



(A.28) 



Again, the factor M 2e is absorbed in the measure for 
convenience. 

We also need to evaluate the various sums over Mat- 
subara frequencies. They can be calculated by a standard 
contour trick, where one rewrites the s um as a c ontou r 
integral in the complex energy plane ijKaoustal Il989|h 
Specifically, we need the following sums 



E 



log[w 2 +u 2 ] = /3w + 21og[l-e- /Jw ] , (A.29) 
1 



E 



LUt + LU 2 



2uj 



[l + 2n(w)] , 



(A.30) 



where n(w) — l/(e l3uJ — 1) is the Bose-Einstein distribu- 
tion function. 

We also need to expand some sum-integrals about 
zero temperature. The phonon part of the spectrum then 
dominates the temperature-dependent part of the sum- 
integral. We can therefore approximate the Bogoliubov 
dispersion relation e(p) by P\/2jl and this gives the lead- 
ing temperature correction. The specific sum-integrals 
needed are 



log [uol + e 2 (p)] 
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